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ABSTRACT 

Observations of deuterated species have long proven essential to probe properties and thermal history of 
various astrophysical environments. We present an elaborated chemical model that includes tens of thousands 
of reactions with multi-deuterated species, both gas-phase and surface, in which the most recent information 
on deuterium chemistry is implemented. A detailed study of the chemical evolution under wide range of tem- 
peratures and densities typical of cold molecular cores, warm protostellar envelopes, and hot cores/corinos is 
performed. We consider two cases of initial abundances, with 1) mainly atomic composition and all deuterium 
locked in HD, and 2) molecular abundances accumulated at 1 Myr of the evolution of a cold prestellar core. 
We indicate deuterated species that are particularly sensitive to temperature gradients and initial chemical com- 
position. Many multiply-deuterated species produced at 10 K by exothermic ion-molecule chemistry retain 
large abundances even when temperature rises above 100 K, and can only be destroyed by dissociation. Our 
model successfully explains observed D/H ratios of many single, double, and triple-deuterated molecules, in- 
cluding water, methanol, ammonia, and hydrocarbons in a variety of environments (cold cores, hot protostellar 
envelopes and hot cores/corinos). We list the most abundant deuterated species predicted by our model in 
different environments of low- and high-mass star-formation regions, as well as key formation and destruction 
pathways for DCO + , DCN and isotopologues of H 2 0, H 3 + and CH 3 OH. 

Subject headings: astrochemistry - molecular processes - methods: numerical - ISM: clouds, molecules - 
stars: circumstellar matter, protostars 



1. INTRODUCTION 

The life cycle of molecules covers a wide range of environ- 
ments in the evolution of gas and dust, beginning from the 
sparse interstellar medium and eventually evolving into stars, 
planets and possibly life itself. Molecular hydrogen cannot be 
easily observed in the cold interstellar medium, therefore, we 
have to use other atomic and molecular tracers. The study of 
deuterium chemistry has proven to be useful to constrain ion- 
ization fraction, temperature, density, and thermal history of 
the interstellar medium ( ISM) (Geiss & Reeves|1981|[Turner| 
[19901 |Crapsi et al.|2005| >. 

More than 160 molecules have to date been observed in 
the interstellar medium, identified using radio and infrared 
observations of their rotational and rovibrational transitions. 
The detected molecules include singly- and multi-deuterated 
species that have be en observed in molecular clouds (ND ; 
Bacmannetal.|2010) , (DCO + , DNC; |van der Tak etaL[ 2009 ), 
hot cor es/corinos (D2CO , HDCO; |Bergman et al.||201lj 
(HD 2 + ; I Vastel et al |2004| i, (DCOOCH3; |Demyk et at'EoTS] 
Margules et al.||2010|l, warm pro t ostellar envelope s (HDO; 
J0rgen sen~& van Disnoeck"l|2010| |Liu et al.||20li| , (DCO + , 
HDC OHFarise et al.|2009T; protopl anetary disks (HD O; |Qi] 
e TaLl|20tel ICeccare'll et aLl|20Q5) , (DCN, DCO + ; [GuT] 



loteau et al . 2006, Qi et al.|2008|irand comets (HDCO; |Kuan| 
et al.|2008|>, (HDO; |Villanueva et al.|2009l|Gibb et al.|2010) , 
(CH 3 D; |Bonev et al |2009||Gibb et al |2010| >. Many deuter- 
ated species remain however elusive due to their low abun- 
dances and/or weak transitions. With the completion of the 
sensitive Atacama Large Millimeter Array (ALMA) we will 
be able for the first time to detect and probe various astrophys- 
ical environments in weak lines of numerous complex and 
rare-isotope molecules, including multi-deuterated species. 
Yet, in order to understand the observations of these 



molecules, feasible models and state-of-the-art chemical net- 
works have to be utilized. Detailed theoretical studies of 
deuterium chemistry is difficult due to the limited number of 
known reactions rates involving deuterated species, and a vast 
amount of hydrogen-dominated reactions. Previous studies 
either made due with what available data there were, or using 



"educated guesses" for missing reaction rates (e .g., Rodgers 
|& Millar|1996HTurner|200Tl|Aikawa et al.|2003) . Also, often 
deuterated species with only one or two D-atoms have been 
considered. First studies of deuterium chemistry included 
only a limited amount of reactions due to restricted computer 



power, such as the one adopted by Brown et al. ( |1988| l, in- 
cluding 93 species connected by 703 reactions. Results were 
accompanied with large errors, but over time the size of the 
chemical networks have increased significantly, and with it 
the level at which we can understand the complexity of the 
chemical evolution in the interstellar medium. 

The isotopologues of H3 + have been known for a long time 
to play a key role in deuterium fractionation in cold environ- 
ments, together with other associated ions such as DCO + and 
H2DO + . The chemical evolution of H 2 D + and the other iso- 
topologues have been found to be closely linked to CO, act- 
ing as the main destruction path ( |Bacmann|2003| l. This leaves 
room for other species to dominate deuterium fraction ation in 
warmer en vironments, such as CH2D + and C2HD + (Millar 
|et al.|1989) . 

The importance of H3 on the efficiency of the deuterium 
fractionation via gas-phase chemistry has been investigated 
by |Roberts et al.| ( |2003] >. They have found that in regions of 
high density and heavy depletion, which is the case in prestel- 
lar cores, inclusion of isotopologues of H3" 1 " significantly en- 
hanced the fractionation of ionic and neutral species. Other 
studies have been done, also proving the importan c e of the 
isotopologues of H3" 1 ", such as |Roberts et al. ( |2003| l, |Roberts| 
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& Mill ar) ([20061) andjvan der Tak| (|2006 ). 

As Roberts et al. (2003 1 have found, the environmen t plays 
an over important role in the chemical evolution. |Roberts 



|& Millar] ( |2000b| l have studied the chemistry a wide range 
of physical parameters, varying density, temperature, initial 
abundances and freeze-out of molecules. They have found 
that fractionation is strongly affected by temperature, and 
agreed reasonably well with the limited data on sulfur-bearing 
species observed in dark clouds. If freeze-out is present, 
molecular D/H ratios can become very high, and gas-phase 
chemistry can result in very large abundances of both singly - 
and multi-deuterated molecules. 



Willacy ( |2007| > and |Willacy & WoodslpOOl)) have studied 
the deuterium chemistry in the outer and inner regions of pro- 
toplanetary disks, and found that in the outer regions, pho- 
todesorption is essential in the midplane for high abundances 
of molecules such as HDO, and that the DCO + fractionation 
ratio is very sensitive to what desorption mechanism is im- 
plemented. In the inner regions, they found their calculated 
D/H ratios on ices to be too high compared to observations of 
comets, which suggests that the model is incomplete, possi- 
bly lacking essential surface reactio ns. Their results are quit e 
different compared to older study of |Aikawa & Herbst] ( |1999| >. 
Aikawa & Herbst have found that their chemical model repro- 
duces higher than cosmic D/H ratios in disks, as observed in 
comets, and that the D/H ratio is sensitive to the adopted ther- 
mal and ionization disk struct ure, regulating D-fr actionation 
via ion-molecule chemistry. Bay et et"aL] ((2010) have con- 
ducted a survey of deuterated species in extragalactic star- 
forming regions, studying the influence of density, tempera- 
ture, far-UV radiation field, cosmic -ray ionization and metal- 
licity on D/H ratios of about 20 deuterated species. Their 
results agree well with the limited observations in external 
galaxies, and provide a list of key deuterated species whose 
abundances are high enough to possibly be detected by Her- 
schel and ALMA. 

The surface chemistry on grains have recently been imple- 
mented into chemical models, as their importance has been 



unveiled as crucial catalysts for complex organics ( Herbst & 
van Dishoeck 2009), and other molecules such as H2 and wa- 
ter. Although having a profound effect on the evolution of 
molecules, it remains poorly understood. Extensive labora- 
tory studies have investigated the importance of surfaces in 
deuterium chemistry (e.g., Hidaka et al. 2006||Na gaoka et al. 
2006), as well as theoretical studies (e.g., iW atanabe |2005[ 
Cazaux et al.|2008| |Kalvans & Shmeld|201 1[ >. As hydrogen- 



and deuterium-bearing molecules are destroyed at approxi- 
mately the same rates, Rodgers & Millar (1996) have found 
that initial high D-fractionation in ice mantles can persist for 
over 10 4 years, which holds true for a wide range of physical 
conditions. This means that it is safe to infer fractionation on 
grains from observations of deuterated molecules in hot cores. 

The Hersc hel observatory is currently gatherin g measure- 
ments (e.g |Vastel et al.|[20T0l |Comito et [20101 |Bac^1 
mann et al. 2010]!, and with the upcoming ALMA, SKA, and 
JWST facilities, the interest in deuterium chemistry will only 
grow. In the advent of the high-sensitivity, high-resolution 
measurements these facilities will contribute, we have con- 
ducted an extensive study of deuterium chemistry in both cold 
and hot molecular environments shielded from ionizing ra- 
diation by the A v extinction of 10 mag. The main aim of 
the present study is to investigate qualitatively and quanti- 
tatively deuterium fractionation and chemical evolution for 



an extended set of diagnostic species, ranging from simple 
mono-deuterated molecules to triple-deuterated complex or- 
ganics, using a newly-developed up-to-date chemical network 
(available online). Another motivation is to try to reproduce 
simultaneously observed D/H ratios of a variety of mono-, 
doubly-, and triply-deuterated species in distinct astrophysi- 
cal environments (e.g., HDO/DCO+/DCN, ND2/D2O/D9CO, 
and CD 3 OH). 

In Section [2] we present the chemical model of Semenov 
|et al. ( |2010| l adjusted to the ISM conditions and use it to 
model chemical evolution in a wide temperature-density pa- 
rameter space. The parameter space represents evolutionary 
stages from cold, dark, low-mass cl ouds to the warm dense 
medium in hot cores/corinos (Section|2.1|i. We give a detailed 



description of the deuterium cloning process and our choice 
of mass-dependent fractionation reaction rates an d br anching 
ratios of relevant dissociation reactions in Section l2~3l In Sec- 
tion[3]this model is used to calculate chemical abundances and 
D/H ratios for assorted species. The general trends as well 
as de taile d e volut ion of assorted species are analyzed in Sec- 
tions 3.1 and|3.2| We compare our results with recent obser- 



vations and theoretical studies in Discussion. Summary and 
Conclusions follow. 

2. MODEL 

2.1. Physical model 

Deuterated species have been observed extensively in va- 
rieties of astronomical objects located locally in the Galaxy 
and at high-z redshifts. In this work we are primarily con- 
cerned with modeling deuterium chemistry in cold low-mass 
prestellar cores and warm/hot high/low-mass cores and cori- 
nos. The evolution of the interstellar medium begins from 
fragmentation of turbulent, partly molecular, partly atomic 
clouds, having temperatures of about 80 K and low densities 
of ~ 10 cm" 3 . Eventually some of the clumps form prestellar 
cores having low temper atures of 8-15 K and moderate den- 
sities ~ 10 4 - 10 6 cm" 3 (|Launhardt et al.||2010[ |Andre et al.| 
[20091 |Snow & McCaII]|2006| l, which is favorable for accu- 
mulating high degrees of deuterium fractionation enabled by 
the ion-molecule chemistry via H3 isotopologues. Eventually 
some clouds begin contracting due to gravitation, increasing 
internal densities and temperatures. Then a protostar is born, 
which begins to heat up its surroundings. The envelope ma- 
terial can t hen reach temperatu res > 100 K with densities < 
10 10 cm" 3 ( |van Dishoeck 2009). The protostellar environment 
then evolves into a protoplanetary disk, which is heavily ir- 
radiated by the stellar FUV and X-ray radiation, and later a 
planetary system. 

In this paper we will concentrate on the evolutionary stages 
ranging from a cold molecular cloud to a hot core, discarding 
the later evolutionary stages. We intentionally choose a wide 
parameter space covering temperatures between 5 - 150 K and 
densities of 10 1 - 10 10 cm" 3 , assuming fixed A v = 10 mag (for 
the standard ISM dust). That is, our study is focused on heav- 
ily obscured regions of the ISM and dark protoplanetary disk 
midplanes. The simulations with the fixed extinction allow us 
to plot the respective abundances and D/H ratios as 2D plots, 
which are easier to interpret and discern general trends. 

2.2. Chemical model 

We have utilized the gas -grain chemical m odel "AL- 
CHEMIC" developed by |Semenov efaT] pOTUl l, where de- 
tailed description of the code and performance is avail- 
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able. The code is optimized for time-dependent evolution of 
large chemical networks including both gas-phase and surface 
species. A brief summary is given below. 

The chemical network is based on the osu.2009 ratefile with 
the recent updates to reaction rates as of Nov 201(Q Several 
tens of pho toreaction rates a re updated using the new calcula- 
tions of van Dishoeck et al. (2006), which are publicly avail- 
ably The self-shielding of H 2 from photodissociation is cal - 
culated by Equation (37) from van Dishoeck & Blake ( 1998 ). 
The shielding of CO by dust grains, H2, and its self-shielding 
is calculated using the precomputed table of |Lee et"aL] ( 1996, 
Table 11). 

We consider cosmic rays (CRP) as the only external ion- 
izing source (as well as CRP-induced FUV photons), using 
the standard CRP ionization rate, Ccr = 1-3 x 10~ 17 s" 1 . The 
gas-grain interactions include sticking of neutral species and 
electrons to uniformly-sized 0. 1 /im dust grains with a stick- 
ing coefficient of 1, release of ices by thermal, CRP-, and 
UV-induced desorption, dissociative recombination and radia- 
tive neutralization of ions on charged grains, and grain re- 
charging. We do not allow H2 and its isotopologues stick- 
ing to grains. Chemisorption of surface molecules is also not 
consid ered. We assume the U V photodesorption yield of 10~ 3 
(e.g., Oberg et al.|2009a" b c >. To allow synthesis of complex 
molecules, an extended list of su rface reactions and photodis- 
sociation of ices is adopted from Garrod & Herbst|( |2006| ). 

We assume that each 0. 1 fim spherical olivine grain provides 

2 x 10 6 sites for surface recombination that proceeds solely 
through the Langmuir-Hinshelwood mechanism. Upon a sur- 
face recombination, there is a 5% chance for the products to 
leave the grain due to release of energy. Following exper- 
imental studies on the fo rmation of molecular hydrogen on 
dust grains by |Katz et aT] ( |1999[ ), we adopt the standard rate 
equation approach to the surface chemistry without quantum- 
mechanical tunneling either through the potential walls of the 
surface sites or through the reaction barriers. 

A typical run, with relative and absolute accuracies of 
10~ 5 and 10~ 15 , the full gas-grain network without deuterium 
chemistry (~ 7000 reactions, ~ 650 species) takes a few sec- 
onds for 5 Myr of evolution (Xeon 2.8GHz CPU). The same 
model with added deuterium chemistry (~ 52,000 reactions 
and ~ 1 , 700 species) takes approximately an order of magni- 
tude longer to calculate. 



2.3. Deuterium chemistry 

As we aim to conduct an extensive deuterium chemical net- 
work, we have fol l owed the same algorithm as described in 
Rodgers & Millar (1996). Namely, reactions bearing hydro- 
gen atoms are considered to have deuterated analogues, and 
"cloned" accordingly (assuming the same rate coefficient if 
no lab data are available). In cases where the position of the 
deuterium atom is ambiguous, we apply branching ratios, di- 
viding the total rate coefficient over the number of reaction 
channels. 

In Equation (H]i a typical example of the cloning process is 
presented, which results in three separate reactions depend- 
ing on where the deuterium is placed. The original reaction 
is to the left, and to the right are the three possible isotopic 
pathways, each with a branching ratios of 1/3. 



C 4 H^ + S^HC 4 S + + H 



C 4 HD + 
C 4 HD + 
C 4 DJ 



+ S 

+ s 
+ s 



DC 4 S + +H 
HC 4 S + + D 
DC 4 S + +D 



See: 


http : //kida . obs . u-bordeauxl . f r 




htt 


: / /www . strw . leidenuniv . nl/~ewi 


ne/photo/ 



(1) 

A few restrictions are applied to the cloning routine. Due to 
the significantly large number of hydrocarbons in the network, 
species with more than six chained carbon atoms (C6-) are 
not isotopically cloned. Any -OH endgroups are conserved as 
well. Observations of deuterated species suggest that fraction- 
ation of species with -OD endgroups is less p robable and can 
be excluded. For example, Pa rise et al.| p006) have conducted 
a survey of deuterated formaldehyde and methanol, and found 
that the fractionation level of CH3OD is at least an order of 
magnitude lower than that of CH2DOH. A hypothesis of rapid 
conversion of CH3OD into CH3OH in the gas-phase due to 
protonation reactions that would affect only species for which 
deuterium is b ound to the electronegative oxygen has been 
suggested by Osamura et al. ( 2004[ >; [Charnley et al. (1997), 
which further support this. 

This method is applied to the latest osu.2009 chemical net- 
work, excluding the recently added fluorine reactio ns, and 
including updated rate coeff i cients fro m Roberts & Millar 
(2000b); Robert s"l!^ ( |2T^ ; |Wakelamet al |p010| ); |Chabot' 
|et al.| ( |2010[ ) as well as from the Kinetic Database for As- 
trochemistrjpl Isotope exchange reactions involving all iso- 
topologues of H t as well as CH2D + , C2HD + a nd others have 
been ta ken from [Roberts & Millaf| ( |2000b| ) and |Roberts et al.| 
d2004| . 

A tabl e of add ed and updated reactions are listed in Ta- 
bles IA1I and IA2I This has resulted in a chemical network 
consisting of 52,548 reactions connected by 1,716 species, 
to our knowledge the most extended network for deuterium 
chemistry to date. 

2.4. Initial abundances 

Reaction rates and physical properties need to be specified 
for the model, as does initial abundances. We have chosen 
to implement two different initial abundances and investigate 
how they affect the chemical evolution over 1 Myr under con- 
sidered parameter space. In both models, all deut erium is ini- 
tially located in HP , with D/H = 1.5 x 10" 5 ( |Stancil et aL 
1998| Linsky|20 03 ). This is valid as the zero-point energies 
of HD and H2 differ by ~ 410 K, h ence favoring production 
of HD over H 2 ( |Roberts et al.|2003| ) . 

For the first model, hence on referred to as the "Primor- 
dial" m odel, we utilized "low metals" abundances o f|Lee et aT] 
(1998), supplied with the HD abundance of 1.5 x 10" 3 (Ta- 
ble[T|). The initial abundances for the second model, the "Evo- 
lution" model, are calculated with our deuterium chemistry 
model using a TMC-1 environment: T = 10 K and «e = 10 4 
cm" 3 , for 1 Myrs. Under such conditions deuterium fractiona- 
tion becomes very effective, leading to high fractionation lev- 
els for many species. The final TMC-1 abundances at 1 Myr 
are used as initial abundances for the "Evolution" model. The 
most abundant deuterated species are given in Table [2] 

3. RESULTS 

We study the general trends of the fractionation ratios of 
deuterated species, looking at dependencies on temperature, 

3 http://kida.obs.u-bordeauxl.fr/ as of [2010-12-20] 
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Table 1 

Initial abundances for the Primordial model. 



Species 


H 2 


H 


HD 


He 


C 


N 


O 




0.499 


2.00 x 10~ 3 


1.50 x 10~ 5 


9.75 x 10~ 2 


7.86 x 10~ 5 


2.47 x 10~ 5 


1.80 x 10 -4 




S 


Si 


Na 


Mg 


Fe 


P 


CI 




9.14 xlO -8 


9.74 x 10 -9 


2.25 x 10~ 9 


1.09 x 10~ 8 


2.74 x 10" 9 


2.16 x 10-'° 


1.00 x 10~ 9 



Table 2 

Initial abundances for the Evolution model. 



Species 


H 2 


H 


HD 


He 


C 


N 


O 




0.500 


1.54 x ur 4 


5.26 x 10" 6 


9.75 x 10~ 2 


6.05 x 10~ 9 


9.18 x 10~ 8 


8.55 x 10~ 7 


Species 


S 

3.13 x 10~ 9 


Si 

1.05 x 10" 10 


Na 

8.64 x l(r" 


Mg 
4.10 x lO" 10 


Fe 

6.36 x 10"" 


P 

1.09 x 10-" 


CI 

1.43 x 10" 10 


Species 


D 

4.20 x 10~ 6 


D 2 

1.86 x 10 -6 


HDO (ice) 
8.64 x 10~ 7 


CH 3 D (ice) 
3.54 x 10~ 7 


OD 

1.27 x 10~ 7 


NH 2 D (ice) 
1.26 x 10~ 7 


CH 3 D 
5.01 x 10~ 8 



density and initial abundance. We also do in-depth stud- 
ies of several assorted species representing certain groups of 
molecules including ions, complex molecules as well as im- 
portant molecular probes. 

3.1. General trends 

Plots of the D/H fractionation ratios for both chemical mod- 
els and assorted species are shown in Figure [T] The left panel 
in each subplot is the Primordial model and the Evolution 
model is shown in the right panel. Species in the four topmost 
subplots show mainly dependence on the initial abundance, 
with differences in fractionation ratios between the models at 
T > 40-50 K. The middle four plots mainly show dependence 
on temperature, and no dependence on the initial abundance. 
The last four subplots show D/H ratios that are sensitive to 
neither the kinetic temperature or the initial abundances. As- 
sorted molecules distributed in these three distinct groups are 
listed in Table [3] 

The computed D/H fractionation ratios reach levels of 
~ 10~ 3 and higher for most of species, which is an enhance- 
ment by ~ 2 orders of magnitude relative to the cosmic D/H 
(~ 1.5 x 10" 5 ). Only at T > 100 K do the ratios reach the 
low cosmic values for some of the molecules. Observation- 
ally this overabundance, primarily in cold environments, has 
been found for many deuterated species. For example, the 
observed abundances are > 10% of their main isotopologues 
for species such as C H 2 DO H, D 2 CO flCeccarelli|2002[ i, D 2 
(|Butner et al.||2007| |, H 2 D + (|CaseUi et al.||2M3E "HBO ( |Liu 
|etal.|2011) and NH 2 D ( |Hatchell|2003l 7eto. We discuss the 
com paris on of observations with our models further in Sec- 
tion |4.1| So far no unique solution has been found to ex- 
plain these deuterium enhancements, though usually either 
gas-phase or grain-s urface fra ctionation routes hav e been in- 
voked (e.g., |Roberts et al.|2004]|Charnley et al.|1997| >. 

The deuterated species are affected differently by variations 
in temperature and the initial abundances. While kinetic tem- 
perature is a key factor for mass-dependent fractionation pro- 
cesses, the choice of initial abundances is also important to re- 
tain certain level of pristine fractionation accumulated at low 
temperatures even when T > 50- 100 K. The sensitivity of 
the D/H ratios to density is, in general, weak since most of 
the chemical reaction rates for main and minor isotopologues 



of a chemical species will be similar. The only exception is 
light-weight multi-deuterated species like D 2 , H 2 D + , D 2 H + , 
etc., having the largest mass differences with their major iso- 
topologues, and thus the largest deviation in their desorption 
and surface reaction rates. It is, perhaps, most evident for 
H 2 D + , shown in Figure [T| and the other minor isotopologues 
of H3 + . These species experience a drop in fractionation ra- 
tios toward lower densities, however only visible towards low 
temperatures. For species that are strongly dependent on the 
production of the H3" 1 " isotopologues, the same dependence 
is noticeable, however somewhat diluted by a larger number 
of intermediate reactions. This can be understood as a mani- 
festation of the unique role of mass-fractionation relevant for 
deuterium chemistry, which is solely dependent on gas tem- 
perature and barely on its density. Therefore we do not treat 
density dependence separately but rather discuss it together 
with the other parameters. 

The majority of species show strong dependency on the ini- 
tial abundances, somewhat fewer strictly to temperature and 
even less, such as CH3D, CHDCO and HDCS, show no de- 
pendency both either parameter. Freeze-out species barely 
show any difference in the distribution of fractionation ratios 
compared to the gas-phase molecules, so generally the same 
general trends apply for species on grains. 

In the "Primordial" model, for species dependent on the ini- 
tial abundance, the D/H fractionation ratios can drop several 
orders of magnitude at higher temperatures (T > 50 - 80 K), 
down to the cosmic value, while in the "Evolution" model 
they usually remain almost as high as in cold environments, 
with only local drops and enhancements. Most radicals are 
strongly sensitive to the choice of the initial abundances, and 
only a few of them show a strict sensitivity to temperature, 
such as DCN and CD 2 CO. Another group of species system- 
atically showing strong dependence on the initial abundances 
are multi-deuterated species. Ions, however, show stronger 
dependency on the gas temperature, as their chemical evolu- 
tion is closely related to that of isotopologues of H3 . Only a 
small number of ions are dependent on the initial abundances, 
such as DCO + and N 2 D + . 

3.1.1. Sensitivity of D/H to the initial abundances 
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Figure 1. Plots of the D/H ratios at 1 Myr for the "Primordial" (left) and "Evolution" (right) models showing different dependencies to physical parameters and 
the initial abundances. The species are (left to right): NH 2 D, H 2 DO+, DCOOH, HDO, D 2 0, H 2 D + , N 2 D + , CH 2 D+, CHDCO, C 2 HD, HDCS, CH 3 D. 
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Table 3 

Species showing dependence on initial abundance, temperature or independent to both. 



Initial abundance 


Temperature dependence 


Independent 


dependence 


Low 


High 




C 2 D, C 2 HD, C3HD, C 4 D, C 4 HD, CD3OH 


C 2 D 2 , CD 3 + 


C 2 D 2 +, C 2 HD+ 


C3H3D, CH 3 D, CHDCO 


CD 4 , CD, CH 2 D 2 , CH 2 DOH, CHD 2 OH, CHD 3 


CHD 2 +, D 3 + 


CD 2 CO, CH 2 D + 


D 2 , CD (ice), HDCS, HD 


D 2 CS, D 2 0, D 2 S, DC 3 N, DC 5 N, DCO + 


D 3 0+, DCN 


H 2 D + 




DCOOCH3, DCOOD, DCOOH, DNC, DNCO, DNO 


HD 2 + 






DOC\ H 2 DO+, HD 2 0+, HDCO, HDO, HDS, N 2 D+ 









ND 3 , ND, NH 2 D, NHD 2 , OD 
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The initial abundances are bound to have a great effect on 
many species, as the abundances of molecules initially present 
in the "Evolution" model allow for deuterated and more com- 
plex molecules to begin forming much faster, and to survive 
for much longer at elevated temperatures. In Figure [1] the top 
four species show strong dependence on the initial abundance, 
as the two models show significantly different distributions of 
fractionation ratios, in particular, at high temperatures. For 
many species, the "Evolution" model has close to a uniform 
distribution of fractionation ratios of > 10~ 3 , while "Primor- 
dial" models show a significant drop toward higher tempera- 
tures (> 100 K). 

Species are affected differently by changes in the physical 
parameters due to the pace at which they are produced, de- 
stroyed, accreted, and evaporated. This in turn is determined 
by the number of reaction pathways, their reaction rates, and 
how many intermediate molecules their formation is depen- 
dent on. In the "Primordial" model, there are initially only the 
basic elements present, along with H2 and HD, which, as soon 
as the simulations begin, are redistributed to form new species 
of increasing chemical complexity. The adopted low T = 10 K 
in this model results in efficient deuterium fractionation via 
ion-molecule chemistry driven by the D-isotopologues of H3 
initially, and other channels later (neutral-neutral gas-phase 
and surface reactions). The reverse processes in this case typ- 
ically imply substantial reaction barriers (> 500- 1 000 K) to 
reach the cosmic D/H abundance values. Consequently, the 
"Evolution" model starts with complex molecules, including a 
large fraction of overabundant deuterated species, that may re- 
tain their abundances even at T > 100 K. Therefore, "Primor- 
dial" model show big differences compared to the "Evolution" 
model, especially, at T > 50 K. A vast majority of the species 
dependent to the initial abundance are radicals, spreading be- 
tween combinations of H, C, O and S-bearing species, up to 
complex organics. A small number of key ions are also found 
in this group: DCO + , DOC + , N 2 D + , H 2 DO + and HD 2 + . 

In Figure |2j we show a few example of multi-deuterated 
species sensitive to the initial abundance. Multi-deuterated 
species are formed in a long chain of formation and D- 
fractionation channels, causing the production speed to be 
slow. In the "Primordial" model there are none of the preced- 
ing species present initially, and producing multi-deuterated 
species becomes time-consuming. In the "Evolution" model, 
in contrast, the preceding deuterated species can be present 
in large amounts to start forming even more complex multi- 
deuterated species early on during the evolution. Therefore, 
many of the multi-deuterated species have different distribu- 
tions of the D-fractionation ratios between the two models. 

3.1.2. Sensitivity of the D/H ratios to temperature 

A number of species shows dependency on temperature, 
without any significant differences between the two models 
of initial abundances (see Figure [T] middle four plots). At 
temperatures exceeding about 20 — 50 K, the fractionation ra- 
tios drop, either abruptly, or smoothly, and both models show 
a similar gradient with temperature. These dependencies trace 
back to either one or several dominant chemical channels be- 
ing activated at certain temperatures, decreasing the abun- 
dance of the deuterated species or slowing down the produc- 
tion. 

Ions represent a majority of the species dependent strictly 
on temperature, although the most common ions are depen- 
dent on the initial abundances. The best example of such 
species is H 2 D + , shown in Figure [T] and the other isotopo- 
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Figure 2. Plots of the D/H ratios at 1 Myr for the "Primordial" (left) and 
"Evolution" (right) models for assorted multi-deuterated species showing de- 
pendence on the initial abundances. The species are (left to right): D?CS, 
CH 2 D 2 and CHD 2 OH. 

logues of H3" 1 ". It shows a strict barrier with a temperature of 
20 K, after which the fractionation ratios drop smoothly and 
reach levels comparable to the cosmic level (~ 10~ 5 ) at about 
100 K. We present further disc ussion of the evolution of the 
isotopologues of H3" 1 " in Section 3.2. 1 

This group of species with strict dependency on tempera- 
ture can be split into two subgroups, diverging by the crit- 
ical temperature where the D/H ratios begin to decrease. 
The deuterated species formed via low-temperature fraction- 
ation channels involving isotopologues of H3 have a criti- 
cal temperature of ~ 40 K, whereas other D-species syn- 
thesized via high-temperature fractionation channels involv- 
ing CH2 D + and C 2 H D + show higher critical temperature of 
~ 80 K ( Parise et al.|2009| l. Several multi-deuterated species 
belong exclusively to the low-temperature group. 

3.2. Detailed chemical analysis for assorted species 

We have looked more thoroughly at a selection of species, 
chosen by their importance as representatives of different 
groups of molecules. Among the selected species, we 
have representative ions (isotopologues of H3 + , DCO + ), sim- 
ple molecules (DCN), complex organics (isotopologues of 
CH3OH) and water (HDO, D 2 0), which is gaining a lot of 
interest now with Herschel observations and the upcoming 
ALMA facility. 

We analyze the evolution of these representative species in 
four distinct astrophysical environments: densities of 10 4 and 
10 8 cm" 3 , and temperatures of 10 and 80 K for the "Primor- 
dial" and "Evolution" models. Using our robust "Chemical 
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Analyzer" (CHAN) tool, we isolate the most important forma- 
tion and destruction pathways for each of these species. Path- 
ways are considered dominant if they contribute to at least 
10% of the net abundance change. Reactions with branching 
ratios are considered important when their combined rate con- 
tributes to > 10% of the net abundance change. We list key 
formation and destruction path ways for the assorted spe cies 
and their main reactants in Table [BT] [B2] [B3] [B4] [B5] and |B6j 
and discuss their formation briefly below. 

3.2.1. H\ isotopologues 

The isotopologues of H3" 1 " have been a target of extensive 
investigation, most recently by, among others, Roberts & Mil- 
lion l|2006), |Sipila et all pOlO) ; |VasteI| (|2007|> ; IVastel et al. 
( 2006b[); |vari der Tak| ( |200^1 i7arid |RoEerts et al.| Q2003| l. There 
is a large interest in these isotopologues as they serve as excel- 
lent probes of physical properties of the interstellar medium, 
and as the starting point for ion-molecule chemistry. The for- 
mation of the H3" 1 " isotopologues is believed to be well un- 
derstood, formed through ongoing deuteration of H3" 1 " and its 
isotopologues by reacting with HD. H3 + is formed through 
H 2 + , which is formed by cosmic rays ionizing molecular hy- 
drogen. Key f orm ation and destruction pathways are summa- 
rized in Table |Bl| and for essential reactants in Table |B2| but 
essentially formation proceeds as follows: 



ion 



H 2 +hi/ = 

H^ + HD: 

H 2 D + + HD = 
HD 2 + HD = 



m 

H 2 D + + H 2 
HD^ + H 2 
D3 + H 2 



(2) 
(3) 
(4) 
(5) 
(6) 



The isotopologues have been observed extensively; e.g. 



H 2 D + (Caselli et al. 2008 ; Cernicharo et al. 2007 1 Harju et al. 



2006] l and HD 2 + ( |Parise et al.|201 l||Vastel et al.|2004| i. While 
observations of the ortho- or para-states of these isotopo- 
logues have been successful, we have not considered it in 
our models. The fully deuterated D3" 1 " is symmetric and not 
observable in the sub-mm range, and the same goes for H3" 1 ", 
making it hard to determine fractionation ratios for the species 
directly. H3 + can instead be measured using IR/optical ab- 
sorption of its ele ctronic transition toward warm environments 
( |Goto et al.|2008] l. 

All the H3 isotopologues show dependence on temperature, 
and partly on density, as can be seen in Figure [3] In general, 
the fractionation ratios reach values above 10 , and only at 
temperatures above ~ 30 K does it go lower. The decrease 
is noticeable as a clear fractionation boundary towards the 
higher temperatures, which is more pronounced for the multi- 
deuterated isotopes. At lower temperatures there is weak de- 
pendency on density, with higher fractionation ratios achieved 
toward higher densities. 

D3" 1 " show significantly higher D/H fractionation ratios at 
low temperatures, increasing with density, up to values as 
high as 10-20, while other H3 isotopologues reach values of 
10 _I - 1 at most. At T < 20 K, backward de-fractionation re- 
actions are not active, and fractionation goes fast because the 
main destruction channel via the ion-molecule reaction with 
CO is absent due to all CO being locked in grain icy mantles 
( Vastel et al.|2006b{ Bacmann et al.|2003~ i. Since the energy 
barriers of the backward reactions cannot be overcome, and 
D3" 1 " cannot be further deuterated, it will be a final product of 
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Figure 3. Distribution of the D/H fractionation ratios over the investigated 
parameter space, for the three isotopologues of H3 + and the two sets of the 
initial abundances. 

the H3 fractionation, and thus it reaches the highest D/H frac- 
tionation compared to those of HD} and H 2 D + . 
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Figure 4. Relative abundances at 1 Myr of the CO molecule, both in gas- 
phase and ices for the "Primordial" (left panel) and "Evolution" (right panel) 
models. 

In the limit where CO is depleted onto grains at T < 20 K 
(see Figure [4}, D3 + can theoretically reach values ~ 20, while 
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H2" 1 " and HD2" 1 " reach levels of unity (Roberts 



, which is 



in agreement with our results. Finally, the drop of the D/H ra- 
tios for the H3 isotopologues observed in Figure|3]at densities 
exceeding 10 8 cm -3 is due to steady removal of the elemental 
deuterium in the gas-phase in deuterated ices, such as HDO, 
D2O, CH3D, NH2D, etc., by surface processes. This reduces 
the pace of the deuterium fractionation via H3 as abundances 
of HD and other gas-phase D-species go down, whereas abun- 
dances of H3 increase, thus decreasing the corresponding D/H 
ratios. Long timescales associated with these surface pro- 
cesses require high densities to have profound effect on the 
H3 abundances during the 1 Myr evolutionary time span. 

An area of high D/H at ~ 20 K, towards densities above 
~ 10 6 cm" 3 , is present for all the H3 + isotopologues, which we 
will refer to as fractionation "islands". Two big nearby peaks 
are evident that reach D/H close to or above unity, especially 
evident for D3 + . As we discussed above, this is due to ex- 
treme molecular depletion of major destruction agent for the 
H3 isotopologues, namely the CO gas, and the evolution be- 
comes restricted to a limited set of reactions. Under these cir- 
cumstances, the network becomes bistable to slight changes, 
because a certain species may be more/less depleted than oth- 
ers. Therefore, it is merely a feature of our model, and since 
the isotopologue radicals are highly reactive, acting as start- 
ing point for the formation of many molecule, it will affect the 
distribution of other deuterated species as well. 

3.2.2. H2O isotopologues 

Water is one of the most essential molecule in the context 
of planet habitability and the evolution of complex lifeforms, 
and is gaining a lot of interest in the light of results on ex- 
oplanets from Kepler, as well as the Herschel results and in 
light of upcoming ALMA facility. A lot of effort is put on un- 
derstanding the presence of water throughout the evolution of 
the interstellar medium and protoplanetary disks. In this con- 
text, we perform a deeper analysis of the chemical evolution 
of the deuterated isotopologues of water; HDO and D2O, and 
deduce their main reaction pathways. 

Plots of the fractionation ratios for deuterated water, both 
gaseous and solid, are shown in Figure [5] The highest mod- 
eled D/H ratio for HDO and D 2 is > 10%, achieved at tem- 
peratures below about 20 K and moderate densities of 10 6 - 
10 8 cm" 3 . There are some significant differences between the 
two water isotopologues in the two initial abundances mod- 
els. The "Primordial" model shows strong temperature de- 
pendence for the both isotopologues, although not as clear 
as for the isotopologues of H3" 1 ". None of the most dominant 
formation reactions for either HDO or D2O show any tem- 
perature dependence, suggesting instead that it is inherited 
from the c hemi cal evolution and fractionation efficiency of 
H3" 1 ". Table B3 lists the key formation and destruction path- 
ways associated with the water isotopologues, mainly formed 
by dissociative recombination of H2DOVHD2O" 1 ". Essentially 
the formation of the water isotopologues proceeds as follows: 



H 2 DO + + e" 
HD 2 + +e" 
CHDCO + +e" 
HD 2 DO + +e" 
D 3 + +e" 
CD 2 CO + + e" 



HDO+H 
HDO + D 
HDO+C2 
D 2 + H 
D 2 + D 
D2O + C2 



(7) 
(8) 
(9) 
(10) 
(11) 
(12) 
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Figure 5. Distribution of the fractionation ratios for the isotopologues of 
H 2 0. 

The listed reactions dominate throughout the considered 
parameter space, but toward higher densities reactions with 
HCN/HNC also become significant. The main reactants 
H2DO" 1 " and HD2O" 1 " are dependent on the isotopologues of 
H3" 1 " through the intermediate formation of HDO + . Thus, the 
temperature-dependence becomes somewhat diluted by the 
intermediate reactions pathways. The decrease in D/H for 
water occurs at a slightly lower temperature than for H2D + , 
due to the longer, slower evolution of HDO and D2O. The full 
list of domina nt re action channels for H2DO + and HD2O" 1 " are 
listed in Table |B2| There are two dominant surface reactions, 
but we note that desorption only plays a small role in increas- 
ing the abundance of the H3O" 1 " isotopologues, only at higher 
temperatures does it become significant. 

In the "Evolution" model, HDO has an almost constant 
fractionation ratio over the whole parameter space, both in the 
gas-phase and on grains. The HDO abundances are quickly 
accumulated in the previous cold evolutionary phase used as 
an input for the "Evolution" model. Since the mass differ- 
ence between HDO and H2O is only 5%, both these species 
stick to and desorb from dust grains at similar physical con- 
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ditions. Consequently, their abundance ratio remains nearly 
constant at a wide range of T and bh, particularly for their 
ices. While HDO is able to retain a high fractionation ratio 
(> 10~ 3 ), even at high temperatures in the "Evolution" model, 
D2O only show a slightly increased D/H ratio at temperatures 
< 20 K. It is no surprise that D2O is harder to produce, as most 
multi-deuterated species have lower abundances compared to 
their singly-deuterated analogues. 

The main destru ction pathways for the water isotopologues 
are listed in Table B3 Overall, the strongest channel for re- 



moving deuterated water is depletion onto grains, which be- 
comes dominant after ~ 1 4 — 1 5 years. Reactions with dif- 
ferent ions such as HCO + , H3 + and C + , while individually less 
important, are able together to make a significant contribution 
to removing water. 



3.2.3. DCO + isotopologues 
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Figure 6. Distribution of fractionation ratios for both gas-phase and iced 
DOC+ and DOC+ in both models. 



DCO + is an important molecular ion, which is used as a 
tracer of deuterium fractionation and ionization fraction, and 
has an isotopomer DOC + . We distinguish between these two 
forms in the chemical network, however not in the analysis as 
we find their evolution similar. The D/H fractionation ratios 
over the investigated parameter space are plotted in Figure |6j 
for both DCO + and DOC + . As can be clearly seen, the iso- 
topomers show similar D/H distributions, where the temper- 
ature gradient is less smooth for DCO + compared to DOC + . 
Over the whole parameter space the general fractionation ratio 
remains at a level of ~ 10~ 3 , reaching values as high as > 1 at 
T < 20-30 K and «h > 10 cm" 3 . The same "island" features 
as in the case of H3 are evident here, which can be attributed 
to the fact that DCO + and DOC + are daughter molecules of 
the H3 isotopologues. This is further strengthened by the sim- 
ilarity in the temperature gradient between these deuterated 
ions. 

In Table |B4| the major formation and destruction reactions 
are listed for the both isotopomers, but essentially the forma- 
tion proceeds through the following pathways: 



H 2 D + + CO^DCO + +H 2 
HD^ + CO^>DCO + +HD 
D 3 + CO^DCO + + D 2 
N 2 D + + CO^DCO + + N 2 



(13) 
(14) 

(15) 
(16) 



The key formation process is an ion-molecule reaction of 
CO with different protonated ions. These include the iso- 
topologues of H3" 1 ", HCO + and N 2 + . Another dominant pro- 
duction pathway is a substitute reaction, where HCO + react 
with deuterium, and CH 2 D + reacting with oxygen. These for- 
mation pathways are only dominant in the cold environments 
where CO is severely depleted from the gas, whereas at higher 
temperatures the CO production channel dominates. The last 
pathway is able to change the structure of the isotopomers by 
reacting with molecular hydrogen, and favoring the formation 
of more energetically favorable DCO + over DOC + , causing 
the slight differences for their D/H distributions in the "Pri- 
mordial" and "Evolution" models. 

Key destructive pathways are listed in Table |B4| however 
the most dominant one is the dissociative recombination with 
electrons, splitting the molecular ions into CO and D. Since 
CO is the most common reactant to form DCO + , it will con- 
tribute to the reformation of DCO + . The same is true for the 
less dominant destruction pathways involving C, HCN/HNC 
and H, which may be counteracted by the increased CO abun- 
dances. 

3.2.4. DCN isotopologues 

Observations of HCN have been extensively used as probes 
of dense molecular gas and photon-dominated regions (PDRs) 
(Boger & Sternberg 2005). Both it and its deuterated isotopo- 
logue DCN has an isotopomer, DNC, and the formation path- 
ways are similar between the two. Observati ons of HCN and 
HNC show an estimated ratio of unity ( |Herbst et al.||2000[ ), 
and laboratory studie s of DCN and D NC estimate this ratio to 
be about unity (Hiraoka et al.|[2006] >. Therefore, in the anal- 
ysis we do not discern between the reaction pathways of the 
two structural configurations of DCN. 

In Figure [7J the fractionation ratios for the gas-phase and 
solid DCN and DNC are plotted. The corresponding D/H ra- 
tios reach levels of > 10% and more, until the temperature 
barrier of ~ 80 K is reached, at which it drops down to the 
cosmic level (~ 10 ). Just as for DCO + , the temperature bar- 
rier is not well defined, possibly due to a dependence on the 
H3" 1 " isotopologues diluted by intermediate reactions. There 
are extended features of high D/H visible at high temperatures 
(T ~ 40-50 K) for the "Evolution" model, most prominently 
for DCN. 

The D/H distribution suggest that the chemical evolution 
of the HCN/HNC isotopologues is strongly coupled to that 
of HCO + and Ht. The key formation pathways are listed in 



Table B5 but essentially the formation proceeds through the 
following reaction pathways: 



H 2 D + + HCN/HNC = 
HD 2 +HCN /HNC = 
HDCN + + e" : 
D 2 CN + + e" : 
HCN/HNC + D = 



-HDCN + + H 2 
•D 2 CN + + H 2 
-DCN/DNC+H 
-DCN/DNC+D 
DCN/DNC + H 



(17) 
(18) 
(19) 
(20) 
(21) 
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Figure 7. Comparing results between the two models for the DCN/HCN 
ratio over the investigated parameter space. 

Neither of the formation reactions include any of the H3 
isotopologues as reactants. Instead, this dependence origi- 
nates through other reactants such as HDCN + and D2CN" 1 " in 
reactions 19 and [20] w hich depend on the H3 + isotopologues 
through reactions |1/| and 18 The H3 + isotopologues react 
with HCN/HNC, further increasing the fractionation ratios. 
DNC/DCN can also form through HNC/HCN reacting with a 
deuterium atom in reaction[2T| but is in general only dominant 
during the first 1 4 — 1 5 years of evolution. 

The HCN/HNC isotopomers are mainly destroyed by reac- 
tions with different ions such as H3" 1 ", HCO + and He + . Rela- 
tive to the three main formation pathways however, they are 
small, and the only effect they have is to slightly slow down 
the chemical evolution. Freeze-out is instead the main chan- 
nel for removing DCN and DNC, but it is only dominant at 
higher densities (~ 10 8 cm" 3 ), while at higher temperature 
(~ 100 K) it is counteracted by thermal desorption. 

3.2.5. Complex organics: methanol isotopologues 

Complex organics are suggested to be the starting blocks 
of more advanced molecules essential in the evolution of life, 



such as amino acids, the structural units for proteins. Many 
of them are also excellent probes of the surface chemistry 
in the mid-plane and molecular layer of protoplanetary disks 
( |Semenov]2 010). Since they form primarily in cold environ- 
ments, they are also of interest for investigation in the earlier 
stages of the interstellar medium. 

In Figure [8] the D/H fractionation ratios of the CH3OH iso- 
topologues are plotted for both models. The D/H distribu- 
tion is not as smooth as for other species which we have pre- 
viously discussed, and we see lower fractionation ratios to- 
ward higher temperatures. The ratios can reach well above 
unity for CH2DOH, especially in "Evolution" model, and or- 
ders of magnitude lower for the other two isotopologues. An 
"island"-like variation of the methanol D/H ratios are also no- 
ticeable in Figure [8] especially for CHD2OH. 

Much as for other species, there are significant differences 
between the two models. The D/H ratio for methanol remain 
well above the cosmic level in the "Evolution" toward higher 
temperatures, though in the "Primordial" model this trend is 
less pronounced. CH2DOH is able to retain fractionation lev- 
els above unity in a substantial region of the parameter space 
in the both models. In the "Evolution" model all multi-D 
methanol isotopologues have D/H ratios that are well above 
the cosmic level, whereas in the "Primordial" model this is 
true only at T < 50 K. 

There is a dependence on the physical properties that affects 
the methanol production pathways. Initially in the "Primor- 
dial" model, up to approximately 10 2 - 10 5 years (depending 
on density), reactions [22] -[24] are the primary formation path- 
ways. Later, the other isotopologues become dominant, such 
as CH 3 OHD + , CH 3 OD 2 + and CH 2 DOHD + , also reacting with 
free electrons. At higher temperatures, desorption becomes 
significant, eventually reaching an equilibrium between des- 
orption and freeze-out. 



CH 2 DOH^ + e _ 
CHD 2 OH^ + e _ 



CH 2 DOH + H 
CHD2OH + H 



CD3OH} + e" CD3OH + H 



(22) 
(23) 
(24) 



In the "Evolution" model, the formation pathways are the 
same but less effective. Instead, desorption and freeze-out of 
pristine methanol dominate. At densities > 10 s cm" 3 freeze- 
out is the strongest, and at temperatures above ~ 70 K desorp- 
tion dominates. The CH3OH isotopologues are only formed 
in the colder (< 50 K) environments, and elsewhere the envi- 
ronment reach a state of equilibrium after ~ 100 years. 

Many small destruction pathways are active throughout the 
evolution, their individual contribution to the net abundance 
change is, however, le ss th an 10%. The most active reac- 
tions are listed in Table |B6| showing that there are not many 
reactions competing in the formation and destruction of the 
CH3OH isotopologues. The same reaction pathways are ac- 
tive in the both models, however at somewhat lower rates in 
the "Evolution" model, where also freeze-out is active as well 
as some minor desorption channels. 

The isotopologues of CH3OH2" 1 " are the reactants of the 
dominant formation pathways, and are formed in a chain of 
reactions mainly consisting of molecules like CH+ reacting 
with H2 or H2O. The un-deuterated equivalent chain of pro- 
cesses is presented in Eq. 25 We see that there is a chemical 
link, albeit weak, betweenthe isotopologues of CH3OH and 
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the H3" 1 " isotopologues. 
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Figure 8. Distribution of fractionation ratios for the isotopologues of 
CH3OH, both in gas-phase and on ices. 



H£+C = 
CH + +H 2 = 
CH^+H 2 = 
CH3+H 2 = 



CH + +H 2 
•CHj+H 
CH^+H 
CH 3 OH 2 



(25) 



The lack of smoothness in the plots may stem from the 
fact that the chemistry of complex organics are not very well- 
known yet, and involve many surface reaction pat hways, part 
of which is mis sing in our chemical network (Herbst & van 
|Dishoeck| |2009). This could also be the "feature" of the 
model, where slight delays with the evolution (or evapora- 
tion) of D-species compared to lighter H-species can cause 
strong alterations in respective D/H ratios over a tiny range of 
temperatures and densities, which show up as local peaks and 
drops of D/H in the figures. 

4. DISCUSSION 

4. 1 . Observation 

Many deuterated species have been observed in the ISM, 
protostellar envelopes, and protoplanetary disks in the recent 
years, and we compare their abundances with predictions of 
our extended deuterium network. A list of observed deuter- 



ated species and their D/H fractions are presented in Table C 1 
where we also list the corresponding model values and the 
spatial scale of the observed objects. We have listed the most 
abundant deuterated species for environments of dark clouds 
(Table Rl, infrared dark clouds (IRDC, Table [5]), high-mass 
protostellar objects (HMPOs, Table [6]l and hot cores (Table 
[7]), the full list is available upon request from authors. 

The modeled values are chosen in an interval of the con- 
sidered parameter space that is typical of the observed object. 
For example, dark clouds are the coldest and less dense en- 
vironments, whereas hot cores/corinos are the hottest and the 
densest. Class -I objects (prestellar cores) typically have tem- 
peratures 10 K and densities 10 4 cm" 3 . It is in these objects 
that the majority of deuterated species have been observed and 
detected. Class objects (warm protostellar envelopes) have 
temperatures > 20-50 K and higher densities of ~ 1 4 — 1 6 
cm" 3 . The Class I objects (hot cores/corinos) have tempera- 
tures > 100 K and densities > 10 6 - 10 8 cm" 3 . 

As there are substantial error bars associated with the accu- 
rate estimation of the physical properties of these objects, e.g. 
due to uncertainties in dust emissivity, temperature, poorly 
known dust-to-gas mass ratio, etc., we take our model val- 
ues from a larger range of temperatures and densities for each 
class of the environments. Also, many of the published frac- 
tionation ratios are based on the measurements made in sev- 
eral sources, and their physical properties can differ signifi- 
cantly, as do the measurements. 

The majority of modeled D/H ratios agree well with the 
observed quantities. Examples are observations in TMC-1 of 
DCO + , DCN, HDS, HDCO, NH 2 D ( |vanDishoecketal.|1995) 
that fit well within an order of magnitude with our model val- 
ues. For a number of species, observations have been incon- 
clusive, and only upper limits are available. An example of 
such is H 2 D + , which has been observed extensively previ- 
ously, but due to the complexity of observing H3 + , there is 
only one available measurement of H 2 D + , merely with upper 
limits. These isotopologues could possibly become observ- 
able with the future ALMA facility. 
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Table 4 

Most abundant deuterated species in cold dark clouds. 



Species Model Observation 



HD 


1.45 x 10- 


-5 


D 


5.04 x 10- 


-6 




3.02 x 10- 


-7 


HDO (ice) 


2.42 x 10- 


-6 


CH 3 D (ice) 


7.94 x 10- 


-7 


NH 2 D (ice) 


3.44 x 10- 


-7 


OD 


2.34 x 10- 


-7 


CH,D 


1.51 x 10- 


-7 


HDO 


5.98 x 10- 


-8 


CH 3 D (ice) 


4.62 x 10- 


-8 



10"' -10- 5 [1,2, 3,4] 
1.8- 2.4 x 10~ 5 [5] 



< 2.5 x 10- 3 [6] 
0.07 [7] 



Referenc es. — (1) Lacour et al. 12005 1; (2) Snow & Mc- 
|CaUl(2T)06l; (3)|Balashev et al.H20l'07r( 4) Bacmann "eTaT 
{2t)(0\; (5)|Rogers et al.H2M7) ; (oT fffllen et al.Hlj)74f ; (7)' 
|Panse et al.|j2UU5) 



Table 5 

Most abundant deuterated species in infrared dark clouds 
(IRDC). 



Species Model Observation 



D 


1.01 x 10-5 


HD 


7.38 x 10-6 


HDO (ice) 


3.68 x 10-6 


D 2 


2.24 x 10-6 


CH,D (ice) 


1.06 x 10-6 


NH 2 D (ice) 


6.38 x 10-7 


D 2 (ice) 


3.38 x 10-7 


CH 2 D 2 (ice) 


2.18 x 10-7 


NHD 2 (ice) 


2.06 x 10-7 


ND 3 (ice) 


1.18 x 10-7 



1.8- 2.4 x 10~ 5 [1] 
10~ 7 -10- 5 [2,3,4,5] 



Ref erenc es. — (l)|Rogers et al.M2007h (2)|Lacour et al. 
2005 ; (3) Snow & McCall|12D06fri4) |Balashev et~aT 
2010 ; (5) Bacmarmeni.l^OiO) 



Table 6 

Most abundant deuterated species in high-mass protostellar 
objects. 



Species Model Observation 



HD 


2.46 x 10- 


-5 


HDO (ice) 


2.26 x 10- 


-6 


D 2 


6.98 x 10- 


-7 


CH 3 D 


6.30 x 10- 


-7 


D 


5.28 x 10- 


-7 


NH 2 D (ice) 


2.52 x 10- 


-7 


D 2 (ice) 


6.96 x 10- 


-8 


CHiD 2 


3.78 x 10- 


-8 


HDOt 


1.83 x 10- 


-8 


NHD 2 (ice) 


1.57 x 10- 


-8 



References. — (1) Teixe ira et aL]jl999) 



For H2D" 1 " and the other isotopologues, we believe to have a 
good understanding of the basic chemistry thanks to extensive 
laboratory studies. However in recent years further studies of 
ortho- and para-modes have been conducted as observations 
have revealed these modes in the interstellar medium. As our 
model does not include the modes of the H3 + isotopologues 
and other species, we are unable to fully compare our models 
with observations. 



Also for later evolutionary stages of the interstellar 
medium, in prestellar cores, and hot cores/corinos, we find 
that the models are in good agreement with observations. A 
few species stand out however. For D2CO, HDCO and NHD2, 
our models underestimate the observed values. When com- 
paring our resu lts of these species to the models of Ro berts &] 
Millar (2000a), we only find a good agreement with NHD2. 
We will discuss this discrepancy further in the comparison 
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Table 7 

Most abundant deuterated species in hot cores. 



Species Model Observation 



HD 


2.32 x 10- 


-5 


HDO (ice) 


3.92 x 10- 


-6 


D 2 


1.05 x 10- 


-6 


CH,D 


4.36 x 10- 


-7 


D 


1.17 x 10- 


-7 


HDO 


1.06 x 10- 


-7 


D 2 (ice) 


5.58 x 10- 


-8 


CH 2 D 2 


2.66 x 10- 


-8 


NH 2 D (ice) 


9.56 x 10- 


-9 


CHD, 


4.24 x 10- 


-9 



0.001-0.01 [1] 



2.94 x 10~ 2 [2] 



Ref erences. — (1) Teixeira et al. 1 1999 1; (2) Butner et al. 
12007) 



with models. 

For a large number of essential species, we have not found 
any published fractionation ratios, and the reason to this may 
be the same as for H3 + , namely a symmetric structure, or the 
signal is too weak to observe, or it has simply not yet been 
found. We can however expect a boom of new observations 
to be revealed with ALMA. Overall, our models are able to re- 
produce the observed fractionation ratios, and we await anx- 
iously new observations from ALMA for new comparisons. 

4.2. Models 

A multitude of models have been used over the years since 
the first molecule was observed in interstellar space, evolv- 
ing to implement larger chemical networks, including multi- 
deuterated species and surface chemistry. Also physical treat- 
ments have improved, such as inclusion of vertical and radial 
mixing, evolution of physical properties and more. 

Incorporating a small network consisting of 2445 gas-phase 
reactions, linked by 225 species, Rodgers & Millar (1996) 
investigated the chemical evolution of the deuterium chem- 
istry in hot cores, with T > 100 K and «e — 10 7 cm -3 . They 
looked especially at the evolution of an assortment of species; 
HDCO, CH3D, NH 2 D, HDO, DC 3 N, DNC/DCN, DCO + , C 2 D 
and OD, which we compared to our model for similar temper- 
ature and density range. 

The same method was incorporated to our chemical net- 
work to incorporate deuterated species , how ever the deu- 
terium chemistry of Rod gers & Millar] ( |1996| l incorporated 
merely 100 deuterated species linked by 1,000 reactions, 
while our networks include approximately an order of mag- 
nitude more. While we incorporated a normal branching ratio 
where reaction rates were spread evenly among reaction path- 
ways created through the cloning, they implemented statisti- 
cal branching ratios, meaning that some species may have in- 
creased or decreased abundances due to the changed branch- 
ing ratios. The results of our models agree well within an 
order of magnitude. This comparison however is limited to 
a small set of species, and there may very well be discrepan- 
cies in the evolution of other species, as our understanding of 
many s pecies has evolved since the study of Rodgers & Millar 
( |1996|) was conducted. 

|Roberts & Millar] ( |2000b| l did a similar study to ours, look- 
ing at the chemical evolution over a parameter space in tem- 
perature and density, however over a somewhat tighter param- 
eter space with lower resolution. They implemented a pseudo- 
time-dependent model that incorporated most reactions from 



the UMIST97 database ( |Millar et al.11997) , consisting of al- 
most 300 species linked by over 5,000 reactions. They only 
included simple surface chemistry for the formation of H 2 and 
HD, while we include a more extensive network of reactions. 

We compare the results of our model for a number of 
species, including DCO + , HDCO, DCN/DNC and NH 2 D, 
looking not only at fractionation levels but also the distribu- 
tion of D/H over the parameter space. We first only compare 
values for a TMC-1 environment after 10 5 years, and summa- 
rize values between our models for several essential deuter- 
ated species in Table [8] We note that the agreement is good 
within an order of magnitude for all but one species; NH 2 D, 
where the fractionation ratios between the models diverge by 
a multiple of ~ 15. It is diffi cult to de termine the source of 
this discrepancy, but Roberts & Millar ( 2000b ) implemented 
the UMIST database, while we base our study on the osu.2009 
network. Due to significant difference in rates at low temper- 
atures between these chemical network, we believe it to be the 
cause of this big discrepancy. 

Looking instead at the parameter space study, and the gen- 
eral trends for the same species, we find the biggest discrep- 
ancy for DCO + and DCN/DNC. For other species the mod- 
els agree well within the ratio intervals of our models, and 
we see similar dependencies to temperature and density. For 
DCO + there are significant differences between the models. 
Where our model keeps a steady fractionation ratio up until 
very high temperatures of ~ 100 K, theirs drop rapidly already 
at 20 K. The same goes also for DCN/DNC, with similar dis- 
tributions as to DCO + . This, as well, may well be caused by 
th e two different netw orks implemented, as discussed above. 
d2000bb 



Roberts & Millar 



co nsidered 

deuterated species, but in a connected paper|Roberts & Millar 



only 



singly- 



(2000a) looked also at multi-deuterated species as well as sur- 
face chemistry. They only discuss a limited set of species, but 
for NH 2 D and the isotopologues of water, HDO and D 2 0, we 

find a good agreement. 

NHD 2 show ratios at or below the cosmic level in Ro berts &| 



Millar (2000a), even at the lowest temperatures, while in the 
same temperature range our model retain levels of 10~ 3 - 10 _I 
up to a temperatur e ~20 K. However, our models agree with 
observations, e.g. |Roueff et al.| ( |2000] |2005 >. Once again it 
is difficult to determine the source of this discrepancy, but we 
believe it to be caused by the different chemical networks im- 
plemented as we discussed previously. In the model where 
surface chemistry is turned off, HDCO also show consider- 
able differences to our model values, which further show how 
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Table 8 

Comparison of fractionation ratios for TMC -1 environment (T = 10 K, 
n H = 10 4 cm" 3 ) between our model and|Roberts~& Millar 12000b}. 



Species 


Our model 


|Roberts& Millar l2000a| 


NHtD 


4 x 10~ 3 


8 4 x 10~ 2 


HDCO 


5.1 x 10~ 2 


4.2 x io~ 2 


DCN 


2.0 x 10" 2 


0.9 x 10 -2 


DNC 


1.9 x 10~ 2 


1.5 x 10~ 2 


C 2 D 


1.6 x 10~ 2 


1.1 x 10~ 2 


C 4 D 


1.4 x 10~ 3 


0.4 x 10~ 2 


DCO+ 


7.4 x 10~ 2 


1.9x 10~ 2 


N 2 D + 


7.6 x 10~ 2 


2.5 x 10~ 2 


C 3 HD 


7.4 x 10~ 3 


0.6 x 10~ 2 


C3H3D 


1.0 x 10~ 2 


8.3 x ur 2 


DC3N 


1.1 x 10~ 2 


0.7 x ur 2 


DC 5 N 


2.7 x 10 -3 


2.3 x 10~ 2 


HDCS 


2.4 x 10~ 2 


4.0 x 10~ 2 



important grain surfaces are as catalysts for organic organic 
molecules. 



In the study of Roberts et al. ( 2004| l, two chemical net- 
works; Rate99 and osu2003, were implemented and com- 
pared as to their effect on the results. They used the model to 
explain recent observations of the CO depletion in prestellar 
cores, and its association to D2CO and HDCO fractionation 
ratios, and found reasonably good agreement. Between the 
two chemical networks, the agreement is good, with fraction- 
ation ratios for most deuterated species agreeing well within 
an order of magnitude. The biggest disagreement occur for 
some of the sulfur-bearing species, where the reaction net- 
works diverge the most. 

Comparing fractionation ratios between our models for 
TMC-1 environments, the majority of species agree well, in- 
cluding H 2 D + , N 2 D + , DCO + and HDO But a number of 
species also show big discrepancies between our studies, 
which dominantly include multi-deuterated species such as 
D 2 0, HD 2 + , D 3 + and NHD 2 . As their study is mainly aimed at 
investigating the effects of implementing the different chemi- 
cal networks, Roberts et al. ( 2004 1 incorporated only a sub-set 
of the networks where only species with six or fewer carbon 
atoms in total were added. Because surface chemistry was 
not included as well, many essential reaction pathways are 
missing, which may lead to the big disagreement between our 
models. This show how sensitive the more complex as well as 
multi-deuterated species are to the networks, and the depen- 
dence on large networks. 

In a more recent study, Roueff et al. (2007 ) used steady- 
state chemical modeling to confirm recent observations of 
high fractionation ratios of DCN/HCN i n dense clumps of the 
Orion Bar, observed by Leurini et al. ( 2006|l. They incor - 
porated the Meudon gas-phase network ( Roueff et al.| 2005 ), 
including 218 species with multi-deuterated isotopologues in- 
volving the elements C, N, O and S, connected by more than 
3250 reactions, but only for gas-phase chemistry. They fo- 
cused on studying the physical conditions of warm clumps, 
and assume them to be homogeneous, with physical param- 
eter ranging in temperature of 10-70 K and densities «h = 
10 4 - 10 7 cm" 3 . They incorporate two different initial abun- 
dance as well, both atomic, where some material frozen onto 
grains in cold source is assumed to be desorbed into the gas in 
one, and the second similar to our adopted cosmic abundance. 

Among the discussed species in their study, the majority 



of species agree well with our models. Three species how- 
ever show greater discrepancies; DCO + , HDCO and D 2 CO. 
For DCO + , there are big differences in the frac tionation ra- 
tios be tween the two incorporated models by Roueff et al. 
( |2007[ ). The discrepancy are dependent on the initial abun- 
dance adopted, and illustrate how species such as DCO + are 
affected by the initial abundances, even to slight changes. 

HDCO and D 2 CO show similar fractionation ratios in both 
model, but the discrepancy stems instead from the lack of sur- 
face chemistry. In our models, at low temperatures, significant 
amounts of HDCO and D 2 CO are formed on grain surfaces, 
and desorbed into the gas phase, being the main source of for- 
mation for the species later in the evolution. 

4.3. Conclusions 

We have shown how our model simultaneously can explain 
observed D/H ratios of both cold and hot environments of the 
interstellar medium, implementing a simple two-staged evo- 
lution model and, to our knowledge, the biggest chemical net- 
work for deuterium chemistry to date. Our model is not with- 
out its drawbacks however, and we discuss these here and the 
effects that they may infer on our results. 

In the past decade there has been a multitude of studies and 
experiments that investigate the different structural c onfigura- 
tions that species may exist in; ortho, para and meta ( Crabtree 



l et al.1201 1 1 |Vastel et al.|20T0l|Pagani et al.|2009||C"aselIi et al. 
|2008||Vastel et al.|2006a|b| l~These modes have been found to 
play a crucial role in the behavior of several species, includ- 
ing H 2 , water, and the isotopologues of H3" 1 ", and if not taken 
into accou nt, it can have significant effects on the chemical 
evolution (Flower et al.|2006| . Given the complexity and size 
of the chemical network, we chose not to include the different 
modes of deuterated and non-deuterated species. 

Our simple evolutionary model, where we simulate a TMC- 
1 environment for 1 Myr and use the final abundances as the 
initial abundance for the model, is a very simple approach to 
the evolution of the interstellar medium during the chemical 
evolution. It is only a 2-staged model, whereas a number of 
previous studies have incorporated a continuous evolution of 
the environment during their evolution. The advantage to our 
approach is that we are able to look at a much wider parameter 
space with higher resolution, however, the simple approach 
of evolutionary model will affect the chemistry through the 
constant change physical parameters. As we have shown in 
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Section [3] the chemical evolution of a majority of species is 
dependent on the temperature, while density only has minor 
effects on a limited number of species. The increase in tem- 
perature may affect the evolution of some species by activat- 
ing new destructive pathways or releasing important reactants 
into the gas-phase, enhancing other formation pathways. 

Less significant effects on the evolution would come 
from our choice of simulating a set of homogeneous zero- 
dimensional objects with varying densities and temperatures. 
For cold, dark clouds this is valid, as they are more or less 
homogeneous, but for hot cores/corinos this is not a good 
approximation anymore, as temperatures especially change 
drastically with distance from the very centers of the objects. 
But our parameter space study has the advantage that we can 
show of the chemistry change throughout the object's struc- 
ture. Would we simulate hot core/corinos specifically, with 
a physical size and temperature and density gradients, we 
would need to include other processes such as mixing in our 
models, which affect the evolution as well. 

Our choice of using normal branching ratios, rather than 
statistical, may cause over- or under-estimations of reaction 
rates for certain pathways. However, since the reaction rates 
of many reaction pathways are uncertain, this effect will not 
be significant, and rates and branching ratios for the most es- 
sential species have already been measured in laboratories. 

These simplifications in our models affect our results, but 
comparing our models with observed D/H ratios we find good 
agreements for many interstellar molecules, including wa- 
ter, methanol, ammonia and many hydrocarbons. We also 
have very good agr e ement with previous model studies by 
URoueff et al.||2007) |Roberts et al.||2004| |Roberts & Millar| 
2000a|b[ |Rodgers & Millar||1996| l, where the most essential 
species agree within an order of magnitude. The biggest dis- 
crepancies all through are for the species DCO + , HDCO and 
D2CO, where the latter two have been better understood after 
the introduction of surface chemistry into chemical models. 
DCO + has been studied extensively both in laboratories and 
with telescopes, and our model agree well with the observed 
values. 

We also list the most dominating formation and destruc- 
tion pathways for DCO + , DCN and isotopologues of water, 
H3" 1 " and CH3OH, and note that especially surface reactions 
are in dire need of updated reactions rates, especially for com- 
plex organics such as the CH3OH isotopologues. With Her- 
schel gathering new measurements from observations, and in 
the advent of the upcoming ALMA and possibly JWST fa- 
cilities, we have listed predicted relative abundances of the 
most abundant deuterated species in low- and high-mass star 
formation environments of the interstellar medium. 

The research leading to these results has received funding 
from the (European Community's) Seventh Framework Pro- 
gramme [FP7/2007-2013] under grant agreement no. 238258. 
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APPENDIX 

A. UPDATED AND ADDED REACTIONS TO CHEMICAL NETWORK 



Table Al Added and updated non-deuterium reactions. 
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13 


7 




CH 2 +H 


=> 


CH + H 2 


2.20 x 10~ 1U 


0.00 
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CN + O 


=S> 


CO + N 


2.60 x 10~ u 


-0.12 





1 


CN + N 


=S> 


C + N 2 


1.00 x 10" 10 


0.00 





1 


NH + O 


=> 


NO + H 


6.60 x 10-" 


0.00 





1 


c 2 +o 


=> 


CO + C 


2.00 x 10" 10 


-0.12 





1 


C2H + 


=> 


CO + CH 


1.00 x 10~ 10 


0.00 





1 


C3H + 


=> 


C 2 H + CO 


1.00 x 10~ 10 


0.00 





1 


NO + N 


=> 


N 2 +0 


4.76 x 10"" 


-0.35 





1 


NH2+O 


=> 


HNO + H 


6.39 x 10-" 


-0.10 





1 


NH 2 +0 


=> 


NH + OH 


7.10 x 10~ 12 


0.10 





1 


HNO + O 


=S> 


NO + OH 


3.77 x 10-" 


-0.76 





1 


NH3+ CN 


=> 


HCN + NH 2 


2.77 x 10 _n 


-0.85 
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C3N + O 


=> 


C 2 N + CO 


1.00 x 10~ 10 


0.00 
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C 2 H + N 


=> 


CN + CO 


1.00 x 10" 10 


0.00 
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C++H 2 


=> 


CH 2 + 


2.00 x lO" 16 


-1.30 


-23 
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C 2 H 2 + + H 2 


=> 


C 2 H 4 + 


2.90 x 10~ 14 


-1.50 





1 


CH 3 ++ H 2 


=> 


CH 5 + 


3.97 x lO" 16 


-2.30 


22 
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H 3 + + O 


=> 


OH++ H 2 


7.98 x lO" 10 


-0.156 


-1.41 
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H 3 + + O 


=> 


H 2 + + H 


3.42 x 10~ 10 
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H 2 CO++ e+ 


=> 


HCO + H 


1.50 x ur 7 


-0.70 
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H 2 CO + + e+ 


=> 


CO + H + H 


2.50 x Kr 7 


-0.70 
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CNC++ e+ 


=> 


CN + C 


3.80 x 10- 7 


-0.60 
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HCNH++ e+ 


=> 


HCN + H 


9.62 x 10~ 8 


-0.65 
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HCNH++ e+ 


=> 


HNC + H 


9.62 x 10~ 8 


-0.65 
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HCNH++ e + 


=> 


CN + H + H 


9.06 x 10~ 8 


-0.65 
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HC 5 NH++ e+ 


=> 


C5N + H 2 


8.00 x 10~ 8 


-0.70 
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HC 5 NH++ e+ 


=> 


HC5N + H 


9.20 x 10~ 7 


-0.70 
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HC 5 NH++ e+ 
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NH 2 +0 
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O + NH 


=> 


OH + N 
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C 2 H + CRPHOT 
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C 2 +H 


4.05 x 10 3 
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C 3 + CRPHOT 
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C 3 +H 


3.25 x 10 3 
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c 3 +c 


7.70 x 10 2 


0.00 





2 


C 3 H 2 + CRPHOT 
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C 3 H + H 
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c 4 + c 
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1.66 x 10~ 7 


-0.50 





2 
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C 2 + CH 2 
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C 6 ++e+ 
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6.00 x 10- 8 
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Q + +e+ 


=> 


c 6 + c 2 


C 9 ++ e+ 


=> 


c 7 + c 2 


C 10 + +e + 


=> 


c, + c 


C 10 ++ e+ 


=> 


c 8 + c 2 


C 2 H + PHOTON 


=> 


C 2 +H 


C 3 H + PHOTON 


=> 


C 3 +H 


C 4 + PHOTON 


=> 


C 3 +H 


C 4 + PHOTON 
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C2+C2 


C3H2+ PHOTON 


=> 


C3H + H 


C3H2+ PHOTON 
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C 3 +H 2 


C 4 H + PHOTON 


=> 


C 4 + H 


C 4 H + PHOTON 
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C 2 H + C 2 


C 5 + PHOTON 
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c 3 +c 2 


C 6 + PHOTON 
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c 5 +c 


C 7 + PHOTON 


=> 


c 6 + c 


C 8 + PHOTON 
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c 7 + c 


C 2 H + CRPHOT 
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C + CH 


C3H + CRPHOT 
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C2H + C 
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C 2 +CH 


C 4 + CRPHOT 
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C 3 +H 2 


C3H2+ CRPHOT 
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=S> 


C 2 H + C 2 


C 5 + CRPHOT 


=> 


C2+C2 


C 6 + CRPHOT 


=> 
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c 3 +c 3 


C 7 + CRPHOT 
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c 4 + c 3 
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=> 


C6 + C2 


C 8 + CRPHOT 
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C5 + C3 


C 8 + CRPHOT 
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C4 + C4 


C 9 + CRPHOT 
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c 6 + c 3 
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C 3 H++ e+ 
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C 3 H 2 + + e+ 
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c 5 + c 4 


Cio + +e + 


=> 
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C5 + C5 
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1.16 x 10"" 9 


0.00 





2 
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9.00 x Iff" 11 


0.00 





2 


1.00 x 10~ n 


0.00 





2 
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9.50 x 10 2 


0.00 





2 
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2.00 x 10 2 


0.00 





2 
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0.00 





2 
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0.00 
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0.00 
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0.00 
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0.00 
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1.00 x 10" 11 
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2 


6.00 x I0' u 


0.00 





2 


6.60 x 10" 10 


0.00 





2 


2.80 x Kr 10 


0.00 





2 


1.14 x 10-" 


0.00 





2 


1.14 x lO" 11 


0.00 





2 


7.98 x Kr 10 


0.00 





2 


3.42 x 10"" 


0.00 





2 


2.50 x 10" 10 


0.00 





2 
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Reaction 








fl 


■*v 
/ 




C 9 + PHOTON 




C 8 + C 





0.00 





2 


C 9 ++ e+ 




c 8 + c 





0.00 





2 


C 9 + CRPHOT 


=*• 


c 8 + c 





0.00 





2 


C 4 H++ e+ 


=> 


C3+CH 





0.00 





2 


H 3 ++ e+ 


=> 


H 2 +H 


1.36 x icr s 


-0.50 





3 


H 3 ++ e+ 


=*• 


H + H + H 


5.44 x 1(T 8 


-0.50 





3 


C 2 H 2 + + H 2 




C 2 H 4 + 


2.30 x 1(T 14 


-1.50 





3 


CH 3 ++ H 2 


=> 


CH 3 OH 2 + 


5.50 x 10 -12 


-1.70 





3 


CN + N 


=> 


N 2 +C 


1.00 x 10 -10 


0.18 





4 


CH 4 + CH 


=*• 


C 2 H 4 + H 


1.06 x 10~ 10 


-1.04 





4 


H 2 +CH 




CH 2 +H 


1.20 x 10~ 9 


0.00 





4 


CH 3 + + H 2 




CH 5 + 


3.78 x 10~ 16 


-2.30 





4 


C 5 H 2 N++ e+ 




HCN + C 4 H 


8.00 x 10 -8 


-0.70 





4 


C 5 N 2 N++ e+ 




HC5N + H 


9.20 x 10~ 7 


-0.70 





4 


C 5 H 2 N++ e+ 




HCN + C 4 H 


4.40 x 10 -7 


-0.70 





4 


C 5 H 2 N++ e+ 




HNC + C 4 H 


4.40 x 10" 7 


-0.70 





4 


C 5 C 2 N++ e + 




HC 3 N + C 2 H 


1.20 x 10~ 7 


-0.70 





4 



References. — (l) |Wakelam et al.| ( [20Tol >; (2) |Chabot et aT1 ( |20Tol ); (3) |Roberts et aT1 (|2004); (4) KIDA database 



Deuterium fractionation in the ISM 



Table A2 Added and updated deuterium reactions. 



Reaction 






a. 


P 


7 


Refs. 


H 3 + + HD 


=> 


H 2 D++ H 2 


1.70 x HT 9 


0.00 


0.00 


1 


H 2 D++ H 2 


=> 


H 3 ++ HD 


1.70 x KT 5 


0.00 


220 


1 


H 3 + + D 2 


=> 


H 2 D++ HD 


3.50 x 10~ 10 


0.00 





2 


H 2 D++ HD 


=> 


H 3 ++ D 2 


3.50 x 10~ 10 


0.00 


63 


2 


H 3 ++ D 2 


=> 


HD 2 ++ H 2 


1.10 x KT 9 


0.00 





2 


HD 2 ++ H 2 


=> 


H 3 ++ D 2 


1.10 x 10~ 9 


0.00 


251 


2 


H 2 D++ HD 


=> 


HD 2 ++ H 2 


8.10 x 10- 10 


0.00 





2 


HD 2 ++ H 2 


=> 


H 2 D++ HD 


8.10 x 10" 10 


0.00 


187 


2 


H 2 D++ D 2 


=> 


HD 2 ++ HD 


7.00 x 10- 10 


0.00 





2 


HD 2 ++ HD 


=> 


H 2 D++ D 2 


7.00 x 10- 10 


0.00 


107 


2 


H 2 D++ D 2 


=> 


D 3 ++ H 2 


7.00 x 10~ 10 


0.00 





2 


D 3 ++ H 2 


=> 


H 2 D++ D 2 


7.00 x 10~ 10 


0.00 


341 


2 


HD 2 ++ HD 


=> 


D 3 ++ H 2 


6.40 x 10- 10 


0.00 





2 


D 3 ++ H 2 


=> 


HD 2 ++ HD 


6.40 x 10~ 10 


0.00 


234 


2 


HD 2 ++ D 2 


=S> 


D 3 ++ HD 


8.70 x 10~ 10 


0.00 





2 


D 3 ++ HD 


=S> 


HD 2 ++ D 2 


8.70 x 10~ 10 


0.00 


159 


2 


CH 2 D++ HD 


=> 


CHD 2 ++ H 2 


1.60 x 10" 9 


0.00 





3 


CHD 2 ++ H 2 


=> 


CH 2 D++ HD 


1.60 x 10"" 9 


0.00 


370 


3 


CHD 2 + + HD 


=> 


CD 3 ++ H 2 


1.50 x 10 -9 


0.00 





3 


CD 3 ++ H 2 


=> 


CHD 2 ++ HD 


1.50 x 10~ 9 


0.00 


370 


3 


CH 3 + + D 2 


=> 


CH 2 D++ HD 


4.40 x 10~ 10 


0.00 





3 


CH 2 D + + HD 


=> 


CH 3 ++ D 2 


4.40 x 10~ 10 


0.00 


400 


3 


CH 3 + + D 2 


=> 


CHD 2 ++ H 2 


6.60 x 10~ 10 


0.00 





3 


CHD 2 ++ H 2 


=S> 


CH 3 ++ D 2 


6.60 x 10~ 10 


0.00 


400 


3 


CH 2 D++ D 2 


=> 


CHD 2 ++ HD 


1.20 x 10~ 9 


0.00 





3 


CHD 2 ++ HD 


=> 


CH 2 D++ D 2 


1.20 x 10" 9 


0.00 


400 


3 


H 2 D++ D 


=> 


HD 2 ++ H 


2.00 x 10~ 9 


0.00 





4 


HD 2 ++ H 


=> 


H 2 D++ D 


2.00 x 10~ 9 


0.00 


550 


4 


HD 2 ++ D 


=> 


D 3 ++H 


2.00 x 10~ 9 


0.00 





4 


D 3 ++H 


=> 


HD 2 ++ D 


2.00 x 10~ 9 


0.00 


655 


4 


CH 3 + + HD 


=> 


CH 2 D++ H 2 


1.30 x 10" 9 


0.00 





3 


CH 2 D++ H 2 


=> 


CH 3 + + HD 


8.70 x 10- 10 


0.00 


370 


3 


D++ H 2 


=> 


H++HD 


2.10 x 10- 9 


0.00 





3 


H++HD 


=> 


D++H 2 


1.00 x 10- 9 


0.00 


464 


3 


C 2 H 2 ++ HD 


=S> 


C 2 HD++ H 2 


1.00 x 10~ 9 


0.00 





5 


C 2 HD++ H+ 


=> 


C 2 H 2 + + HD 


2.50 x 10~ 9 


0.00 


550 


5 


D++H 


=> 


H++D 


1.00 x 10~ 9 


0.00 





6 


H++D 


=> 


D++H 


1.00 x 10~ 9 


0.00 


41 


6 


H 3 ++D 


=> 


H 2 D++ H 


1.00 x 10~ 9 


0.00 





7* 


H 2 +H 


=> 


H 3 ++D 


1.00 x 10~ 9 


0.00 


632 


7 


HCO++ D 


=> 


DCO++ H 


1.00 x 10~ 9 


0.00 





7 


DCO++ H 


=> 


HCO++ D 


2.20 x 10- 9 


0.00 


796 


7 


N 2 H++ D 


=> 


N 2 D++ H 


1.00 x 10- 9 


0.00 





7 


N 2 D++ H 


=> 


N 2 H++ D 


2.20 x 10~ 9 


0.00 


550 


7 


OH + D 


=S> 


OD + H 


1.30 x 10~ 10 


0.50 





8 


OD + H 


=> 


OH + D 


1.30 x 10 -10 


0.50 


810 


8 


C 2 H + D 


=> 


C 2 D + H 


5.00 x lO" 11 


0.50 


250 


9 


C 2 D + H 


=> 


C 2 H + D 


5.00 x 10" 11 


0.50 


832 


9 


HCN + D 


=> 


DCN + H 


1.00 x 10~ 10 


0.50 


500 


g* 


DCN + H 


=> 


HCN + D 


1.00 x 10~ 10 


0.50 


500 


9 » 


H 2 D++ e" 


=> 


HD + H 


1.20 x 10" 8 


-0.50 





10 


H 2 D++ e" 


=> 


H 2 +D 


4.20 x 10- 9 


-0.50 





10 


H 2 D++ e - 


=S> 


D + H + H 


4.38 x 10~ 8 


-0.50 





10 


HD 2 ++ e" 


=> 


HD + D 


4.20 x l(r 9 


-0.50 





11 


HD 2 ++ e" 


=> 


D 2 +H 


1.20 x 10~ 8 


-0.50 





11 


HD 2 ++ e - 


=> 


D + D + H 


4.38 x 10~ 8 


-0.50 





11 


D 3 ++ e" 


=> 


D 2 +D 


5.40 x 10" 9 


-0.50 





12 


D 3 ++ e" 


=> 


D + D+D 


2.16 x 10~ 8 


-0.50 





12 


CH 3 ++ D 2 


=> 


CH 3 D 2 + 


3.50 x 10- 14 


-1.00 





12 


CH 2 D++ H 2 


=> 


CH 4 D+ 


2.00 x 10~ 14 


-1.00 





12 


CH 2 D++ HD 


=> 


CH 3 D 2 + 


3.50 x 10~ 14 


-1.00 





12 


CHD 2 ++ H 2 


=> 


CH 3 D 2 + 


3.50 x 10~ 14 


-1.00 





12 


CD 3 + + H 2 


=> 


CH 2 D 3 + 


6.30 x 10~ 14 


-1.00 





12 


C 2 H 2 ++ H 2 


=> 


C 2 H 3 D+ 


3.39 x 10~ 14 


-1.50 





12 


CH 3 ++ HDO 


=> 


CH 3 OHD+ 


1.10 x 10"" 


-1.70 





12 


CH 3 + + D 2 


=> 


CH 3 OD 2 + 


1.65 x 10 -11 


-1.70 





12 


CH 2 D + + H 2 


=S> 


CH 2 DOH 2 + 


1.10 x 10-" 


-1.70 





12 


CH 2 D++ H 2 


=> 


CH 2 DOH 2 + 


1.65 x 10-" 


-1.70 





12 


CH 2 D++ D 2 


=> 


CH 2 DOD 2 + 


2.20 x 10-" 


-1.70 





12 
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Reaction 






a 




P 


7 


Refs. 


CHD 2 + + H 2 




CHD 2 OH 2 + 


1.65x10- 


1 1 


-1.70 





12 


CHDi + + hdo 


=> 


CHD 2 OHD + 


2.20 x 10' 


1 1 


-1.70 





12 


CHD 2 + + D 2 


=4- 


CHD 2 OD 2 + 


2.75 x 10~ 


1 1 


-1.70 





12 


CD 3 + + H 2 


=*• 


CD 3 OH 2 + 


2.20 x 10" 


1 1 


-1.70 





12 


CD 3 + + HDO 




CD 3 OHD+ 


2.75 x 10" 


1 1 


-1.70 





12 


CD 3 + + D 2 




CD 3 OD 2 + 


2.75 x 10" 


1 1 


-1.70 





12 


CH 3 + + HD 




CH 4 D+ 
CH,D 3 + 


0.00 




0.00 





11 


CH,D + + D 2 




0.00 




0.00 





11 


CHD 2 ++ HD 


=*• 


CH 2 D 3 + 


0.00 




0.00 





11 


C 2 H 2 + + H 2 


=> 


C 2 H 3 + 


0.00 




0.00 





11 



References. — (l)|Ceccarelli & Dominik| ( |2005]l; (2)[Giles et al.| ( [T992]l ; (3) |Smith et aL](|1982a|b|l; (4)|Walmsley et al.|(|2004|); 
(5)|Herbst et al.|dl987); (6)|Watson|(|1976|l; (7)|Adams & Smith (1985 ); (8)ICroswell & Dalgarno|(|1985|l; (9)|Schilke et al.l 
( 1992); (lO)^undstrom et al. l994| ; (*! l) |&)berts et al.| ( |2004| ; (lij lOrsson et al.|dl^7T ~ 



Estimate 



Deuterium fractionation in the ISM 



Table Bl 

Most essential formation and destruction pathways for deuterated H3 + ; H2D + , HD2 + and H3 + . 





Reaction 




Rates 


Accuracy 


H 3 ++ HD 
H 3 ++D 
H 2 D++ HD 
H 2 D++ D 
Hi++ D 2 
HD 2 + + HD 
HD 2 + +D 
H 2 D++ D 2 


=*• 
=*• 

=> 


H,D++ Hi 
H 2 D++ H 
HD 2 + + H 2 
HD 2 + + H 
HD 2 + + H, 
D 3 ++ H 2 
D 3 + +H 
D 3 + + H 2 


1.70E-09 O.OOE+00 O.OOE+00 
l.OOE-09 O.OOE+00 O.OOE+00 
8.10E-10 O.OOE+00 O.OOE+00 
2.00E-09 O.OOE+00 O.OOE+00 
1.10E-09 O.OOE+00 O.OOE+00 
6.40E-10 O.OOE+00 O.OOE+00 
2.00E-09 O.OOE+00 O.OOE+00 
7.00E-10 O.OOE+00 O.OOE+00 


(C)(1) 
(C)(2) 
±15% (L) (3) 
(4) 

±15% (L) (3) 
±15% (L) (3) 
(4) 

±15% (L) (3) 


H 2 D++ CO 
H 2 D++ CO 
HD 2 + + CO 
HD 2 + + CO 
D 3 ++ CO 
D 3 ++ e" 


=> 
=*> 
=> 
=> 
=*• 
=*• 


DCO++ H 2 
HCO++ HD 
DCO++ HD 
HCO++ H 2 
DCO++ D 2 
D + D + D 


3.22E-10 O.OOE+00 O.OOE+00 
3.22E-10 O.OOE+00 O.OOE+00 
3.22E-10 O.OOE+00 O.OOE+00 
3.22E-10 0.00E+00 O.OOE+00 
3.22E-10 O.OOE+00 O.OOE+00 
2.16E-08 -0.50E+00 0.00E+00 


< 25% (M)(5) 

< 25% (M)(5) 

< 25% (M)(5) 

< 25% (M)(5) 

< 25% (M)(5) 

< 25% (M)(5) 



Referenc es. — (l) |Sidhu et aT|(T992); (2) | Adams & Smith (1985); (3) |Giles et aE](T992) ; (4) |Walmsley] 
|eTaTH2004) ; (5) Cloned, original from UDFA:] www . udf a . netl 



Laboratory measurement 
c Calculated 

DOMINANT FORMATION AND DESTRUCTION PATHWAYS FOR DEUTERATED SPECIES 



Alberts son et al. 



Table B2 

Important reactions for species involved in the main pathways of the assorted deuterated species; HDO, 
D 2 0, DCO+, DCN, CH 2 DOH, CHD 2 OH and CD3OH. 



Reaction Rates Accuracy 



H 2 + + H 2 




H 3 ++H 


2.10E-09 


0.00E+00 0.00E+00 


<25%(M)(1) 


H 2 + h^cR 


=> 


H 2 + + e- 


9.30E-01 


0.00E+00 0.00E+00 


factor 2 (L)(l) 


HDO++ H 2 


=> 


HoDO + + H 


4.07E-11 


0.00E+00 0.00E+00 


<25%(M)(1) 


D 2 + + H 2 


=*• 


HD 2 + + H 


4.07E-11 


0.00E+00 0.00E+00 


<25% (M)(l) 


H 2 D++ 




HDO++H 


6.68E-11 


-1.56E-01 -1.41E+00 


<50% (M)(l) 


HDo++ O 


=> 


D 2 ++H 


6.68E-11 • 


■1.56E+00-1.41E+00 


< 50% (M)(l) 


H 3 + +C 


=> 


CH++ H 2 


2.00E-09 


0.00E+00 O.OOE+00 


factor 2 (L)(l) 


CH++ Ho 


=*> 


CH 2 ++ H 


1.20E-09 


0.00E+00 O.OOE+00 


<25% (M)(l) 


CH 2 + + H 2 
CH 3 ++ HD 


=*> 


CH 3 + + H 


1.20E-09 


O.OOE+00 O.OOE+00 


<25% (M)(l) 


=> 


CH 2 D++ H 2 


1.30E-09 


0.00E+00 0.00E+00 


(M)(2) 


H,D++ HCN 


=> 


HDCN++ Ht 


1.13E-09 


-0.50E+00 0.00E+00 


<25%(M)(1) 


H 2 D++ HNC 


=> 


HDCN++ H 2 


1.00E-09 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 


HDi++ HCN 


=> 


DoCN++H 2 


1.13E-09 


-0.50E+00 0.00E+00 


<25% (M)(l) 


HD 2 ++ HNC 


=> 


DoCN++ H 2 


1.00E-09 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 


CH 2 +N 


=> 


HCN + H 


3.95E-11 


1.67E-01 O.OOE+00 


<50% (C)(l ) 






uwr 1 


3.95E-11 


1.67E-01 O.OOE+00 




CH + H 2 


=> 


CH 2 +H 


1.20E-09 


0.00E+00 0.00E+00 


< 25% (L)(l) 


CH 3 ++ HDO 


=> 


CH 3 OHD+ 


1.11E-11 


-1.70E+00 0.00E+00 


(3) 


CH 3 + + D 2 




CH 3 OD, + 


1.11E-11 


-1.70E+00 0.00E+00 


(3) 


CH 2 D++ H 2 
CHD 2 + + H 2 


=*• 


CH 2 DOH,+ 
CHD 2 OH,+ 


1.11E-11 
1.65E-11 


-1.70E+00 0.00E+00 
-1.70E+00 0.00E+00 


(3) 
(3) 


CD 3 ++ H 2 


=> 


CD 3 OH 2 + 


2.20E-11 


-1.70E+00 0.00E+00 


(3) 


CH 2 + + H 2 
CH 3 + + H 2 


=> 


CH 3 + + H 


1.20E-09 


0.00E+00 0.00E+00 


< 25% (M)(l) 




CH 5 + 


3.78E-16 


-2.30E+00 2.15E+01 


factor 2 (M)(l) 


CH 3 + + HD 


=*• 


CH 2 D + + H 2 


1.30E-09 


0.00E+00 0.00E+00 


(M)(2) 


CH 2 D + + H 2 
CH 3 ++ D 2 


=*• 


CH 4 D+ 
CHD 2 + + H 2 


2.00E-14 -1.00E+00 0.00E+00 
6.60E-10 0.00E+00 0.00E+00 


(M)(2) 
(M)(2) 


CHD 2 + + H 2 


=*> 


CH 3 D,+ 


3.50E-14 


-1.00E+00 0.00E+00 


(M)(2) 


CD 3 ++ H 2 


=*• 


CH 2 D 3 + 


6.30E-14 


-1.00E+00 0.00E+00 


(M)(2) 


D 3 ++ CH 2 


=^ 


CD 3 ++ H 2 


5.19E-11 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 



References. ( 1) Cloned, original from UDFA: w ww.uc!fa.net| (2) |Smith et"aT]jl982b|a) ; (3) 
|Roberts et al.|(2004) 



M Laboratory measurement 
c Calculated 
L Litterature 



Deuterium fractionation in the ISM 



Table B3 

Most essential formation and destruction pathways for deuterated water; HDO and D2O. 



Reaction Rates Accuracy 



112 i-'w + e 


=> 


MIA' + 1 1 


i one ns 


n ^riFj-fin n nncj_nn 

-U.JUIi+UU U.UUjj,+UU 


< ZJ /C {Wl )\L) 




=> 


nUU + D 


z.ZUc-Uo 


-UJUt+UU U.UUo+UU 


< Zj /C fyL) 


\\i\Jk\) + c 


=> 


nUU + L^ 2 


1 .uuc-u / 


-UJUt+UU U.UUo+UU 


fn / .tnr if \t"\\ 

iaClOr L W^jyY ) 


H 2 DO + + HCN 




HDO + H 2 CN+ 


5.47E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


H 2 DO + + HNC 




HDO + H 2 CN + 


5.47E-10 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 


H 2 DO + + C 3 




HDO + C 3 H+ 


4.00E-10 


0.00E+00 0.00E+00 


factor 2 (L)(l) 


H 2 +OD 


=> 


HDO + H 


9.33E-14 


1.00E+00 1.04E+03 


<25%(M)(1) 


HDO 


=> 


HDO (ice) 








HDO (ice) 


=> 


HDO 








OD (ice) + H (ice) 


=> 


HDO (ice) 








OH (ice) + D (ice) 


=> 


HDO (ice) 








HDO + HCO+ 


=> 


H 2 DO + + CO 


4.20E-10 


-0.50E+00 O.OOE+00 


<50% (M)(l) 


HDO + H 3 + 


=> 


H 2 DO + + H 2 
H 3 0++ HD' 


1.67E-10 


-0.50E+00 O.OOE+00 


<25% (M)(l) 


HDO + H 3 + 




1.67E-10 


-0.50E+00 O.OOE+00 


<25%(M)(1) 


HDO + C+ 




DOC++ H 


6.00E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


HDO + C + 




DCO++ H 


2.97E-10 


-0.50E+00 0.00E+00 


<25%(M)(1) 


HDO + C + 


=> 


HOC++ D 


6.00E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


HDO + C+ 


=> 


HCO++ D 


2.97E-10 


-0.50E+00 O.OOE+00 


<25% (M)(l) 


un,n+ . „- 

n 1 * ' + c 




I \ ( ) , IT 

UlU + 11 


9 OS 
Z.zuc-Uo 


n sriF-i-nn n noFj-nn 

-U.JUETUU U.UUIJ.+UU 


<- Z J /c \ Nl )\LJ 


T-rr*i~f~i+-i_ [ n \ 
rlijiy) + xiciN 




D 2 + H 2 CN+ 


^ a 1 n in 

J. 4 ICj- 1U 


-u. juc+uu u.uuri+uu 


< Z J /c [Jyl )( L ) 


HD20 + + HNC 
HD 2 0++ C 3 




D 2 + HiCN + 


4.95E-10 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 




D 2 + C 3 H+ 


4.00E-10 


O.OOE+00 O.OOE+00 


factor 2 (L)(l) 


D 3 0++ e" 


=*• 


D 2 + D 


2.20E-08 


-0.50E+00 0.00E+00 


<25% (M)(l) 


CD 2 CO++ e~ 


=*• 


D 2 + C 2 


1.00E-07 


-0.50E+00 0.00E+00 


factor 2 (L)(l) 


OD + 2 D 




D 2 + 2 


4.00E-11 


O.OOE+00 O.OOE+00 


< 25% (L)(l) 


OD + OD 


=> 


D 2 + 


5.50E-13 


1.14E+00 5.00E+01 


< 25% (L)(l) 


D 2 


=> 


D 2 (ice) 








D 2 (ice) 


=> 


D 2 








D 2 + C + 


=> 


DOC++ D 


6.00E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


D 2 + C + 




DCO++ D 


2.97E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


D 2 + H 3 + 




HDoO + + H 2 


1.67E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


D 2 + H 3 + 


=*• 


HoDO++ HD 


1.67E-10 


-0.50E+00 0.00E+00 


<25%(M)(1) 


D 2 + H 3 + 


=*• 


H 3 + + D 2 


1.67E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 


D 2 + H+ 




D 2 0++ H 
HDO++ D 


8.11E-10 


-0.50E+00 0.00E+00 


<25%(M)(1) 


D 2 + H+ 




8.11E-10 


-0.50E+00 0.00E+00 


<25% (M)(l) 



References. — (1) Cloned, original from UDFA: www . udf a . net 
M Laboratory measurement 
L Litterature 
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Table B4 

Most essential formation and destruction pathways for DCO + . 



Reaction Rates Accuracy 



H 2 D + + CO 




DCO++ H 2 


3.22E-10 


0.00E+00 0.00E+00 


<25%(M)(1) 


HCO++ D 


=> 


DCO++ H 


1.00E-09 


O.OOE+00 0.00E+00 


(C)(2) 


HD 2 + + CO 


=> 


DCO++ HD 


3.22E-10 


0.00E+00 O.OOE+00 


<25% (M)(l) 


D 3 ++ CO 




DCO++ Do 


3.22E-10 


O.OOE+00 O.OOE+00 


<25%(M)(1) 


N 2 D++ CO 


=> 


DCO++ No 


8.80E-10 


O.OOE+00 0.00E+00 


<25%(M)(1) 


DC07++ CO 


=> 


DCO++ C0 2 


2.47E-10 


-0.50E+00 O.OOE+00 


factor 2 (L)(l) 


DOC++ H 2 


=> 


DCO++ H 2 


1.11E-12 


O.OOE+00 0.00E+00 


<50%(M)(1) 


CHiD++ 


=> 


DCO++ Hi 


4.10E-11 


0.00E+00 0.00E+00 


<25%(M)(1) 


CH 4 D++ CO 
C^HD++ O 


=> 


DCO++ CH 4 


3.51E-11 


-0.50E+00 O.OOE+00 


<25% (M)(l) 


=> 


DCO++ CH 


1.67E-11 


O.OOE+00 0.00E+00 


<25% (M)(l) 


cb + o 




DCO++ e + 


2.00E-11 


0.40E+00 0.00E+00 


< 50% (C)(1) 


DCO++ e" 




CO + D 


2.40E-07 


0.69E+00 0.00E+00 


<25% (M)(l) 


DCO++ SO 




DSO + + CO 


3.30E-09 


-0.50E+00 O.OOE+00 


factor 2 (L)(l) 


DCO++ H 


=> 


HCO++ D 


2.20E-09 


O.OOE+00 7.69E+02 


(C)(2) 


DCO++ C 


=> 


CD++ CO 


1.10E-09 


0.00E+00 0.00E+00 


factor 2 (L)(l) 


DCO++ HCN 


=> 


HDCN++ CO 


2.43E-09 


-0.50E+00 0.00E+00 


<25%(M)(1) 


DCO++ HNC 


=> 


HDCN++ CO 


2.21E-09 


O.OOE+00 0.00E+00 


<25%(M)(1) 



References. — (1) Cloned, original from UDFA: [www . udf a . net) (2) |Adams & Smith|{l985) 
M Laboratory measurement 
c Calculated 
L Litterature 
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Table B5 

Most essential formation and destruction pathways for deuterated DCN. 





Reaction 






Rates 


Accuracy 


HDCN + e 


=> 


DCN + H 


6.17h-08 


-0.65h+00 O.OOh+00 


(1) 


n rw+i ~- 
U2*— ™ ^ ^ 


=r- 


nc v _i_ n 


f. 1 7c no 


n fi^pj-nn n nnFj.no 




riL.JN + JJ 


=> 


i \' , i_i 
DCJN + ri 


1 AAC 1 n 


n eric i nn ^ aae i no 
-U.jUJixUU j.UUiixUz 


(2) 


\TTJ |~\ , p\T 


=> 


DL.IN + INri2 


Z. /O-b-lz 


- 1 . 14tixUU U.UUiixUU 


< JV) /c {Ni)\ j ) 


/~ , urif~ , XT+x_ a— 


=> 


r\/"*KI . 

JJL.JN 4- L.rl 


1 .UUb-U/ 


n ^n/c i nn n nncxnn 
-U.jUrixUU U.UUiixUU 


tactor L j J 


i vn a_ xt 

[J\ \> x IN 




UL^IN + W 


1 vnF 1 n 
1 . /Ut- 1U 


n nnFx.no o oofj_oo 

U.UUiixUU U.UUiixUU 


fn/.t™- 9 CT 

iacior z ^jjj^jj 


v ri 1 J x IN 




|"\f"V , II 
L/L-1N + rl 


1 ITT? 1 1 

1 . jZIi- 1 1 


i ^7c m n nnFj.no 

1.0/ti-Ul U.UUiixUU 


< jU ?o 


rll'v^lN x rl 




FlPT\J -i- Ff„ 
Lyl^lN x rl2 


1 1 1 F 11 


n nnFj.nn n nnFj.nn 

U.UUiixUU U.UUiixUU 


IdCLOl Z UYL^jJ 


v n i 1 'v ,N x C 




nPXT j- PFU 


f. nfYF OR 
O.UUIi-Uo 


n sofj-oo n nnFj-nn 

-U.JUlixUU U.UUiixUU 


IdCLOl Z \\-i)\p) 


C~ T-TP>+_t_ XT 




JJL- IN x V. rl 


o. J Jti- 1Z 


n noFj_oo o oofj_oo 

U.UUiixUU U.UUiixUU 


< Zj 70 (IVlj^J ^ 


r\rtxr+j- ft., 

v^UlN x rl2 




\J{ . \ x v-2'^2 


9 99F 1 ^ 
Z.ZZH- 1 j 


n nnFj.nn n nnFj-nn 
u.uuiixuu u.uutixuu 


<. ZJ /o ULVl^J ) 




=4* 


T^r^XT /icol 
UL-IN ^ICcJ 








DCN (ice) 


=> 


DCN 








DCN + H 3 + 


=> 


HDCN++ H 2 


1.13B-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + H 3 + 




HiCN++ HD 


1.13E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + HCO+ 




HDCN++ CO 


2.43E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + He + 




CN++ He + D 


6.90E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + H + 




DCN++ H 


5.56E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + H 




HCN + D 


1.00E-10 


-0.50E+00 5.00E+02 


(2) 


DCN + H 3 + 




HDCN++ H 2 


5.47E-10 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + H 3 0+ 




H 2 CN++ HDO 


5.47E-10 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + D 3 + 




D,CN++ D 2 


1.13E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + C+ 




CNC++ D 


4.75E-09 


-0.50E+00 O.OOE+OO 


< 25% (M)(3) 


DCN + C + 


=> 


C 2 N++ D 


4.75E-09 


-0.50E+00 O.OOE+OO 


factor 2 (L)(3) 



References. — (1) Cloned, original from KID A: |http : //kida . obs . u-bordeauxl . f r/| (2) 
|Schilke et al.1(T992) ; (3) Cloned, original from UDFA: [^7'. udf a . net| T3) 
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Table B6 

Most essential formation and destruction pathways for deuterated CH3OH; CHiDOH, CHD2OH and 

CD 3 OH. 





Reaction 






Rates 


Accuracy 


CH 2 DOH 2 + + e" 
CH30HD++ e" 
CH30D 2 + + e" 


=> 
=> 
=J> 


CH 2 DOH + H 
CH^DOH + H 
CH 2 DOH + D 


7.38E-10 
7.38E-10 
7.38E-10 


-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 


< 25% (M)(l) 
<25% (M)(l) 

< 25% (M)(l) 


C&DOH 
CH 2 DOH (ice) 


=> 
=> 


C&DOH (ice) 
CH 2 DOH 








CHI>OH 2 + + e~ 
CHtDOHD + + e - 
CH30D 2 ++ e - 


=> 
=> 
=> 


CHD 2 OH + H 
CHD 2 OH + H 
CHD 2 OH + H 


7.38E-10 
7.38E-10 
7.38E-10 


-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 


<25%(M)(1) 

< 25% (M)(l) 

< 25% (M)(l) 


CHD^OH 
CHD 2 OH (ice) 


=> 


CHD 2 OH (ice) 
CHD 2 OH 








CD3OIV+ e" 
CHDtOHD + + e~ 
CH,DOD 2 + + e~ 
CHD 2 OD 2 ++ e" 


=> 
=> 


CD3OH + H 
CD3OH + H 
CD3OH + H 
CD3OH + D 


7.38E-10 
7.38E-10 
7.38E-10 
7.38E-10 


-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 
-0.67E+00 O.OOE+OO 


< 25% (M)(l) 

< 25% (M)(l) 
<25% (M)(l) 

< 25% (M)(l) 


CD3OH 
CD3OH (ice) 


=> 

=>■ 


CD^OH (ice) 
CD3OH 









References. — (1) Cloned, original from UDFA: www . udf a ■ net| 
Laboratory measurement 
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LISTING OF OBSERVED DEUTERATED FRACTIONATION RATES IN COLD, DARK, INTERSTELLAR 

ENVIRONMENTS 



Table CI Listings of observated interstellar deuterated species. 



Species 


Sources 


Spacial scale 


Keis 


Model 




l_ldSS -1 


^IdSS KJ/l 


Beam size [ ] 






U 1 rl 


<^ M-.O A LKJ 




1.0 


i 
1 


1 A~ 3 i a-2 




o 14- n 4 y i n -5 

Z. JIH U.H- A 1U 






z 






1 Q I y( _ 1 A— 5 

1.0 - Z.4 X 1U 






J 






<v. U. 14- 






A 




rlD 1 rl 2 


U. /4 - o.o X 1U 




JUX3U 


c 

J 


< 1U 




6.6 - 70 x 10 






6 






5.98 x 10~ 5 






7 






2.48 - 5.94 x 10~ 5 






7 




ND/NH 




0.3 - 1.0 


41 


8 


io- 4 -io° 


OD/OH 


< 2.5 x lO- 3 




1.8° 


9 


It)" 2 - 10° 


CCD / CCH 


0.01 




33 


10 


l(r 3 -10- 2 




0.045 




33 


10 






0.01-0.18 




20 


11 




D 2 0/H 2 




5 x 10~ 5 


1.5x1.5 


12 


< io~ 5 -io- 3 


D 2 S/H 2 S 


< 0.264 




31 b 


13 


lo-* -ur 2 


/ nv. in 


o nox-0 m 5 






14 


10~ 3 — 10 _1 




0.012 - 0.11 






15 




dco + / Hro + 


013-0 023 




20 


11 


10- 2 -io 




0.02 - 0.18 




20 


16 






0.021 




7x5 


17 






0.006 - 0.04 




25-57 


18 






0.045 ±0.014 




30-96 


19 






0.007 - 0.081 






15 






< 0.03 






4 








0.0006 




14 


lo^-nr 2 






6.1 x 10~4 


10-30 


20 








0.007-0.011 


10-30 


20 




DNC / HNC 


0.015 - 0.03 




20 


1 1 


10" 3 -10 _1 




0.02 - 0.09 




~ 20 


21 






0.008-0.122 




17-20 a 


22 




H 2 D+/ H 3 + 




< 3 x 10~ 3 


13 b 


23 


io^-io -2 


HDO / H 2 


0.07 


>0.01 


10-30 


24 


lo- 3 - ur 1 




< 0.002 


0.03 


10-33 


25 








2.94 x Ur 2 


1.5x1.5 


26 


io~ 2 -ur 1 






< 6 x 10~ 4 


3.1 x2.5 


27 


lO^-lO" 2 


HDO / H 2 (solid) 




0.005 - 0.02 




28 


10" 3 -10 _1 






10- 3 -ur 2 




29 




HDS / HiS 






20 


1 1 


10~ 2 -10 _1 




U.Uo - U.JJ 






It) 






0.03 - 0.04 




26.4 11 


30 






0.017 - 0.052 




~30 


31 






0.04 - 0.44 




11 a 


32 






0.11 




44 


33 






0.03-0.1 




10-20 


34 






0.11±0.03 




18 a 


35 




-- 


~0.1 




9-26 a 


36 








0.042-0.271 


11-16 


37 


10- 4 -ur 1 






0.005 - 0.014 


11-26 a 


38 








0.029 - 0.27 




39 








0.015 


9-26 a 


40 




D 2 CO / H 2 CO 


0.143 ±0.019 




27 a 


41 


io _3 -io- 2 




0.01-0.1 




17 


42 






2.94 x lO" 1 




20-60 


43 






< 0.07 




22 


44 








0.046 - 0.44 


10-30 


45 


io- 5 -io -3 






0.02 - 0.4 




46 








0.01 - 0.04 


20-60 


47 




D 2 CS / H 2 CS 


0.333 






48 


io -3 -ur 1 


HDCO / H 2 CO 


0.014 - 0.015 




20 


11 


io^-kt 1 




1.4 x ur 1 




20-60 


43 






0.005-0.11 




33 


10 






0.107 ± 0.015 




27 a 


49 






0.01 -0.14 




22 


44 








0.094- 1.7 


10-30 


45 


io^-io- 2 






~ 0.15 




50 








0.01 - 0.07 


10-80 


47 








0.14 


20-60 


43 








0.006 




14 





30 
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Species 


Sources 


Spacial scale 


Refs 


Model 




Class -I 


class U/i 


Beam size ['' | 








HDCS / H 2 CS 


0.020± 0.005 




>60 


51 


10 - 


10 




0.333 






48 


10~ 2 - 


ur 1 


ND 3 / NH 3 


~ 8 x 10- 4 




25 


52 


< 10- 5 


-10- 3 




0.024 - 0.033 




22 b 


53 


- 








9.35 x 1CT 4 


25 


54 


< 10- 5 


-10- 4 


NH 2 D / NH 3 


0.02 - 0.1 




20 


16 


10- 3 - 


ur 1 




0.025 - 0.18 




18 


55 


- 






0.10-0.28 




22 b 


54 


10- 4 - 


ur 2 




< 0.02 




20 


11 


- 






0.1 


0.06 


22 


44 








u. 1 - u.o 


/in 
< u. 1 


■7 
/ 


3D 










n n/n s nil 


20 


57 


in-5 
1U — 


1U 






U.U4 - U.JJ 


37 


58 


in-2 
1U — 








0.06 


90 fin 


MO 






r*inLJ2 » 1*1 n 3 


5 x 10~ 3 




22 




10"" 4 - 


ur 1 




fi 05 - 99 
yj.yjj w,__ 




22 b 












0.02 - 0.4 




JU 


< 10 


1 n-3 
-10 


C 3 HD / C 3 H 2 


0.05 - 0.15 




1.7' b 


59 


10 - 


1 n-1 

10 1 


C 4 D / C 4 H 


0.0043 




1.7' b 


60 


10 - 


10 




0.0043-0.018 




17-28 


61 


10 - 


10 z 


C 4 HD / C 4 H2 




0.03 


17-28 


61 


10 - 


10 ' 


CD 3 OH/CH 3 OH 




0.014± 0.014 


15 


62 


< 10 3 


-10 2 


CH 2 DCCH/CH 3 CCH 


0.081 -0.170 




40-60 


63 


10~ 3 - 


ur 1 




0.05 - 0.06 




27 a 


64 


- 






0.07 - 0.10 




40 


65 


- 




CH 2 DCN/CH 3 CN 




> 0.005 


27 a 


66 


io~ 2 - 


ur 1 


CH 2 DOH/CH 3 OH 


0.05 - 0.30 




40-60 


63 


10 1 - 


-10 1 






0.43 - 0.65 


10-30 


45 


10~ 3 - 


10" 1 






0.90 ± 0.3 


11-30 


67 


10" 1 - 


-10° 


-- 




0.37± 0.38 


17-28 


61 


10" 1 - 


-10° 






0.05 - 0.30 




68 


" 




CH 3 (JD/ CH 3 UH 


< 0.1 




20 


1 1 










ni 6 9 s 


in in 
1U-9U 


A ^ 










0.04 ± 0.02 


1 1 7fl 


O / 


lO- 2 - 


-10° 


CHD 2 OH/CH 3 OH 




0.2± 0.1 


11-30 


67 


10- 3 - 


-10° 


DC 3 N/HC 3 N 


0.05-0.1 




38-52 a 


69 


10- 3 - 


10- 1 






0.031 ± 0.011 


17-28 


61 


10" 4 - 


lO" 2 


DC5N/HC5N 


0.006-0.016 




1.4 a 


70 


10- 3 - 


lO" 2 


DCOOCH 3 / HCOOCH 3 


0.02 - 0.06 






71 


10- 3 - 


ur 1 




~ 0.15 




9-33 


72 
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a Half-power beam width (HPBW) 
b Full-width half-maximum (FWHM) 
c LWRS 



